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SYITOPSIS 


This investigation deals with mathematical modelling 
and identification of a batch ball mill grinding system. The 
study has been subdivided into j (1) a pseudo-similarity 
solution to the int egro-diff erential equation of batch 
grinding , (2) exact mathonatical definitions of the parameters 
in the discrete size grinding models in terms of the basic 
size continuous functions , (3) realizations of various 
specialized discrete size models , (4) analysis of some special 
grinding phenomena , (5) identification of the grinding system 
in terms of a mill function , (6) estimation of size discretized 
parameters , and (7) empirical correlations for variation of 
size-discretized rate parameters with particulate mass and ball 
size . A brief summary of the results is given below s 

1. A Pseudo- similarity Solution to the 

Int egro-diff erential Equation of Batch Grinding 

A solution to the int egro-diff erential equation of batch 
grinding in continuous size and time has been obtained using 
Kapur’s similarity solution and standard mathematical transfor- 
mation techniques . The equation describes the couplets 
trajectory of the particle ^ectra in time , exrcept for a small 
initial period , The non- self -preserving nature of the resulting 
size ^ectra has been illustrated by simulation studies. This 
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is ilie most general analytical solution to the grinding equation 
available at present » 

2. Exact Mathematical Definitions of the 
Parameters in the lEscrete Size Models 

Starting from the fundamental int egro-diff erential 
equation of batch grinding , exact mathanatical difinitions of 
the parameters in the two types of discrete size models f cur- 
rently being used , have been obtained in terms of the size 
continuous selection and breakage functions which can be defined 
unambiguously . It hag been shown that , even when the system is 
inherently linear and time invariant, the size discretized 
parameters are , in general, implicit functions of grinding time. 
She time va.riability of the parameters is inherent in the dis- 
cretization procedure itself as the parameters are now functions 
of the particle size distribution within the individual size 
intervals. 

A new function, which plays a central role in the analy- 
sis of discrete size models, has been identified and named 
’Seldcage Function^ L(v,x) = S(v) B(v,x) . Its important pro- 
perties have been delineated . 

3. Realizations of Various Specialized 
Discrete Size Models 

When the particles within a size interval are uniformly 
distributed in size, the most commoiiLy used Bass model is real- 
ized . If the selekage function D(v,x) can be reduced to 12^3?), 
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the ESP model is obtained, which generates the Eo sin-Eammler 

type or Ed sin-Eammler distributions. Given that the sele'fcage 

function is separable i. e. , lCv,x) = I^(v) SSP model 

is obtained. Model equations have been derived and the time 

dependence of the size discretized parameters S- and B. . has 

1 ly 3 

been investigated. The parameters S- and B. . are time -invariant 

1 If 2 

in Bass and RSP models. The SSP model has time dependent and 
time- invariant B. . , Time variation of S. paiameters and the 
inner breakage parameters b. . in SSP model has been computed 
using Kapur's similarity solution. It has been shown that the 
inner- breakage parameter is not ne^igibly small even when the 
sieve size ratio is only ^2 . 

4. Analysis of Some Special Grinding Phenomenal^ 

^ The phenomena of first order disppearance kinetics, 

zero-order production of fines, Alyavdin kinetics, difference 
similar (normalizable) breakage parameters, and variation of the 
S^ parameters with particle size acooj^iing to a power law, have 
been investigated in light of the structure of the specialized 
discrete size models proposed in this study j„ It has been empha- 
sized that a given discrete size model possesses concurrent 
properties none of which must be violated, otherwise serious 
inner- contradict ions will develop. 

5. Identification of the Grinding System 
in Terms of a Mill-Punotion 

A new approach has been proposed for identification of 
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a linear, t ime- invariant , continuous size and discrete time 
grinding systan, in which the size reduction occurs by inpact 
fracture, chipping and abrasion. Batch ball mill grinding data 
for dolomite has been successfully used for illustration of the 
identification procedure . The role of the various grinding 
mechanisms has been discussed in light of the results obtained 
and, the principal advantages of the proposed identification 
scheme have been discussed, 

6, Estimation of Size Discretized Parameters 

Assuming that the size discretized parameters are 
sensibly constant over a narrow time interval, a scheme has been 
devised for exact back-calculation of all tf\e. size discretized 
parameters S- and B. . using Reid’s solution to the batch 

1 1, j 

grinding equation. The type of grinding data required is single 
size feeds ground for a fixed period of time . Next using the 
estimated mill function for dolomite, the required data has been 
generated over several successive intervals of fixed grinding 
time and, for first time, the time variation of the model 
parameters due to their dependence on the size distribution of 
the particles within the individual discrete size intervals, 
has been successfully conputed, 

7, Empirical Obrrelations for Tariation of 
Size~Iliscretized Rate Parameters With 
Particulate Mass and Ball Size 


Sinple systematic gr^hical procedures for determination 



xxx 


of the eii5)irical correlations for the variations in lete 
parameters of the discretized batch grinding equation (the 
Bass model) with particulate mass and the size of the grinding 
media balls have been developed and successfully demonstrated 
for two actual grinding systems . The previously known corre- 
lations for rate parameters have been shoY^n to be restricted 
cases of those obtained in the present study. 



CSaPTSR 1 


miRODUCJTIOfi' 


Cbminiiiution or size redaction of solid materials 
is one of the most important unit operations in many 
metallurgical , chonical and ceramic industries from both 
the technical , as well as , the economic points of view . 

In mineral dressing plants and cement industry , especially, 
each year hundreds of millions of tons of various solids are 
subjected to the size reduction operation , Ihe typical values 
of power consumptions for the required size reduction of 
cement , lead-zinc ore , copper ore , carborundum and 
pignents are 50 , 18 , 12 , 42 , and 66 kWh/ton, respectively 
j la] . In view of very large power consumption in these 
operations , there exists considerable economic incentive for 
improving the design , operation and control of the 
comminution systems . 

In general three broad stages may be identified in 
a size reduction operation , In order of increasing average 
fineness of the feed material these are typically - (i) crushing 
in jaw , gyratory or cone -crushers , (ii) grinding in rod 
mais , and (iii) ball milling , Depending on the initial 
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size of tile feed , its hardness and the size of the product 
desired , in practice all the three steps may iMt be necessary 
in a particular size reduction operation , However ,^in 
many industries the size reduction operation invariably ends 
with ball milling , and whenever it is so , almost the 
entire final product is generated in the ball mills , This 
implies that the ball milling is the most crucial and 
important of all the three steps in the size reduction 
operation , The importance of the ball milling operation is 
further enhanced by the fact that invariably the power 
consumption in the ball mills is several times greater than 
that required for both the preceding size reduction stages 
combined . In view of these ranarks it is not therefore 
surprising that the tumbling mills in general and ball mills 
in particular , have attracted the maximum attention of the 
comminutionists 4 

Till recently , much of the design and operation of 
the ball mills was based on trial and error experience , 
However , ?/ith growing realization of the limitations 
inherent in the empirical approach , there now exists 
considerable interest in the f undemental understanding and 
analysis of the various physical phenomena which operate 
in the ball milling process * e, g. , the rate and mode of 
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breakage of the particles , the size distributions of the 
resulting progeny particles and also the transport of the 

i 

material throu^ the mill , The goal of these efforts is 
that , in principle , it should be possible to fomulate 
a complete quantitative description of the grinding process 
starting from the first principles in terms of only f undaneaital 
■physical p 1X5 perries of the mill - material system and 
various operating variables such as the mill speed * the 

feed rate to the mill and per-cent water in feed etc • The 
process of ball mill grinding being physically quite complex, 
such a development , however , will take considerable more 
effort and time . In order to provide an alternative basis 
for design , analysis , sinulation , optimization and 
control of the grinding process , the conminution engineers 
have in last two decades developed several simplified 
mathonatical models for identification of the mill-material 
system , Depending on the extent to which various physical 
and process variables have been lumped together , these 
models exhibit considerable variations in the approach of 
formulation , the level of sophistication and the details 
of the p2?ocess • Evidently, more sophisticated a model is , 
greater is the number of parameters which characterize 
the model , and therefore , the ta^ of system identification 
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is correspondingly more difficult . On the other hand , over- 
simplified mo dels have very limited utility from the point of 
view of process analysis in the actual practice . 

The present study is based on an approach which is 
now widely employed by the research workers in the field 
of comminution . It is assumed that at any instant the number 
of particles of a given size selected for breakage is 
directly proportional to the total number of particles of that 
size . This is commonly known as the assumption of 'first 
order disappearance kinetics ' . The second assumption is 
that the grinding behaviour of the particles is" not affected 
by the constitution of the particulate asaanbly . In view of 
these two assumptions the syston model in continuous size 
and time variables is therefore , linear in the mass 
fractions of the particles of different size (i. e. , in the 
mass frequency function ) and unless the mill-material 
system or the operating conditions change with time , it is 
also inhere^ntly time invariant , The model incoipo rates 
two basic functions - (i) the selection function , which 
gives the rate of selection for breakage of particles of a 
given size , and (ii) the breakage distribution function , 
which gives the fraction of progeny particles finer than 
a pai^icular size when the parent particle undergo es a single 
primary breakage event « 
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In view of the fact that the experimental size 
distribution data is available only in terns of discrete 
size fractions , size discretized fom of the above model is 
more convenient and also widely used . Hov\/ever , the size 
discretized models currently anployed by various schools 
have either been fomulated without any reference to the 
more fundamental size continuous model and the two basic 
defining functions , or are based on such assumptions v/hich 
are not valid in most practical cases . This arbitrariness 
has led to some controversy and , at the same time , has 
given rise to some interesting problars . First , many such 
discrete size models can be shown to lead to several inner 
contradictions , and not sui^risingly have contributed to 
misinterpretation of numerous experimental data . Second , 
because the parameters of the size discretized model are 
erroneously considered to be completely uncorrelated and 
independent of each other , estimation of a very large 
number of parameters (for example 78 for 12 size intervals) 
becomes extrenely involved in both the experimental as well 
as the computational efforts . 

vVitli these problaas in view , following tasks have 
been attempted in the present study s 
(i) An ai^alytical pseudo -similarity solution to the 
*f undam ental* integro- differential equation of batch grinding 
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has been obtained . 

(ii) Exact mathanatical definitions of the parameters of 
discrete size models in teims of the two basic size continuous 
functions have been established using a fomal mathanatical 

pro cedure , 

(iii) Various specialized discrete size models have been 
identified and the important characteristics of their 
parameters have been investigated . 

(iv) Various special grinding phenomena have been investigated 

in light of the structure of the specialized discrete size 
models , 

(v) A new approach has been proposed and successfully 
danonstrated for identification of the grinding systan in 
continuous size and discrete time . 

(vi) A schane has been proposed for exact back - calculation 
of the parameters using single size feed grinding data , 

The time variation of model parameters due to their dependence 
on the particle size distribution within the discrete size 
intervals , has been computed , 

(vii) G-raphical procedures have been developed and 
danonstrated for determination of the anpirical correlations 
for the variations in rate parameters with particulate mass 
and the size of the grinding media balls . 



CEiPTER 2 
LITEMTURE EEVIM 


In last 50 years or so considerable progress has been 
made in the development of mathanatical models of the size 
reduction processes, Ibr convenience in pres entatioUj, the 
extensive literature on this topic can be classified into four 
broad categories : 

1, Models with a sin^e time- variable parameter 

2, Models based on the concqpt of discrete r^etitive 
breakage steps 

3, Size and time continuous particle population balance 
models 

4, Size discreti 25 ed time continuous particle population 
balance models 

As mentioned in Chapter 1, the discrete size models 
in the last category are of more practical Importance, 
Therefore, a detailed literature review of the parameter 
estimation techniques and empirical correlations for 
variation of these parameters with operating variables has 
also been included. 

2. 1 Models With A Single Time Variable Parameter 

The earlier attempts for the analysis of size 
reduction processes were based on a description of the 
system dynamics in teims of only a sin^e characteristic 
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parameter of tiie particle size distribution as a function of 
the energy received by the particles for breakage j 1-9 1 . 

Typical of such models is the one due to Walker et 
al, |i| who proposed the following relationship 


dx 


£i 

X 


( 2 . 1 ) 


where E is the total specific energy of breakage applied to 
the particles during the course of grinding, r and are 
constants for a given mill-material system, and x is the 
characteristic size of the particle size distribution defined 
as 

r-l- 


1 


Z (M^/x^^') 


; I >~1 


Z M. 
i ^ 


( 2 . 2 ) 


In defining x the entire population of particles has been 
divided into several size groi^js of narrow size interval 
width, each size interval being characterized by an average 
size is the mass fraction of particles in the i-th 

size interval. When 3?=1, the model in Eq. (2,1) is in 
accordance with the Zick*s law of comminution. If r is 
assigned a value 2, the model represents the Eittinger^s 
law, iny other value for r should be regarded only as an 
anpirical constant. 

Due to a lack of proper appreciation of the correct 
definition of x, as given in Eq. (2.2), Walker et al,ts 
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model has been subjected to considerable mi sinteip ret a ti on in 
the subsequent literature. Seme researchers |5-7j have 
confused x for the size of an individual particle and others 
|2-4j 8- 10 1 have assigned to it various arbitrary meanings. 

It should be noted that Bq, (2,2) is not meaningful 
when r=1, As shown by Austin )ll|, in this case correct 
definition of x in continuous size variable should be 


log (x) = f log (x) M (x) dx 


(2.3) 


0 


where M(x) dx is the mass fraction of particulate charge 
in an infinitesimal size range x to x+dx , and Xq is the 
maximum particle size. Similarly the size continuous version 
of Eq, (2.2) is 

Xq 

M(x) dx , r 5 ^ 1 (2.4) 

X ^ ' 0 x^ ' 


or 

X = ^ T ^ i (2.4a) 

where p-. is the (l-r)~th manent of the particle size 
distribution. 

Bond {2,3j noted that the energy irput increases 
as particle size decreases and postulated that the ^eeific 
energy ii^jufe is inversely proportional to a chai^ct eristic 
particle size raised to an exponent 1/2,.- The equation 
proposed by Bond is 
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® “ ^2 C *■ %r5 ~ ”^ 0.5 ] (^* 5 ) 

where subscripts p and f refer to the product and the feed, 
respectively, and is a mill-material constant, This 
principle he has stated as - * ’the energy irput is directly 
propoirtional to the new crack length produced’ * , The implied 
relationship between the characteristic size of the particle 
size distribution x , which he defined arbitrarily as 80 per- 
cent — by weight - passing size, and the crack length, 
however, does not sean to be evident, 

Heplacing the constant 0.5 in Eq. (2.5) by a 
variable parameter r^ , we get the Holmes' model j4l . 

Holmes called r^ as Kick’s law deviation exponent. But in 
view of the difference in the definitions of x in Eqs. 

(2.3) or (2.4) and that assumed by Holmes (i.e. 80 per-cent 

1 

passing size), this inteipretation does not make any sense, 

Charles |5| assumed that the particle size distribu- 
tion of the ground product could be approximated by the 
Gate-Gaudin-Schumann distribution t 

P(x) = ( fr )“ (?»6) 

where E(x) is the fraction of material finer than size x, 
x’ is the size modulus and m'is the distribution modulus. 
Combining Eq, (2.6) with an arbitrary version of Bq, (2. 1) 
in which x is replaced by particle size x, Charles obtained 
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for large reduction ratios 

E= - -f -- ; ; (2.7) 

(3>- 1 ) (m-iM- 1 ) 

It should be pointed out that Bq. (2, l) has no meaning whai 
X is replaced by x jllj , Accordingly, the derivation of the 
Charles' model is untenable. 

Seme workers have formulated their models in terms 
of the ^ecific surface area ^Q* of the particles. Tanaka |8| 
pointed out that with increasing fineness of the material, 
there is progressive depletion of the flaws in the material 
and also the frictional losses increase. Accordin^y, he 
proposed the following anpirical equation 

~ ^4 (Qoo -- Q ) (2.8) 

where Q is the i^ecific surface area , is the limiting 
value of Q, and 0^ is a mill-material constant. 

Harris jsj has preferred to use 'the net energy 
input to the charge* ' E' as the independent variable in his 
model for fine grinding operation in batch mills instead of 
'the energy responsible for fracture inducing eaviiroment' 

E. He argued that as particles become finer dissipetion of 
energy at the points of contact between particles and the 
inter particle pulp increases and, therefore, E is a decreasing 
function of E', He assumed 
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SS = f(E>) = (2.9) 

dE» ^ 


Using the logistic growth hypothesis he proposed the following 
expression 


1 f§ = K - P (Q) 


( 2 . 10 ) 


where K, the potential growth factor, ' according to 
Harris, »‘a constant in conformity with the Kick's law. Ihis 
stat an ent^ however, is incorrect. The expression for Q is 

Q = a / 1 M(x) dx (2.11) 

P 0 

Oomhining Eq, (2,11) with Eg, (-2.3) and substituting in 
Sg. ( 2 . 1 ) withf—l, it can easily be confimed that the 
expression 


Q dE ^ 

does not follow from the Kick's law. 

The second teim E(Q) on EHS in Eg. (2, lO) has been 
temed as 'retardation factor' by Harris, E(Q) is an 
increasing function of Q and has been introduced to account 
for the increasing resistance to fracture as particles 
become smaller , The boundary condition on E(Q) is 

5 (Qj = K (2.12) 

and the assumed functional fonn for F(Q} is 

m = 


( 2 . 13 ) 



Combining Eq.s. (2.9), (2.10) and (2.13) Harris obtained 


dQ KQ( 1 — Q/(^x>) . 

di. — F 


which integrates to 

E 1-(Q/Qoo) 


(2,15) 


where E is the energy to p 2 X)duce half the maximum surface 
area 


It should be pointed out that all the models 

discussed above have a common major drawback. She only 

energy input to the system v/hich can be directly measured 

is the gross energy input E , which is given by 

& 

t 

E„ = / P(t) dt (2.16) 

^ 0 

v/here P(t) is the power draught to the mill. However as the 
models have been formulated either in terns of B or E* , it 
is obvious that the validity of such models cannot be 
confirmed without making further assumptions about the 
relationships of E and E' with , 

Harris [ 12[ has also proposed an alternative model 
now using grinding time as the independent variable, Ihe 
model is descriptive and is based upon a detailed analysis 
of the energy utilization in tumbling mills and various 
grind limit mechanisms , Pormally, the model can be 
r^ reseated by the following equation 
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<i(e/e_ ) 

-dT^ = 2 (e/e^ ) 


(2.17) 


where € is some measure of the fineness (e. g. , specific 
surface area ), is the limiting value of the fineness 
and Z is the product of five different unspecified functions , 
Substituting for Q, and in Eq. (2.14) by 6 , 6^ ani t, 
respectively, and using the following equation 


"m 




(I)" 


(2. 18) 


which, according to Harris, provides best fit to many grinding 
data, the following equation was obtained 


_ s , ^ i 

- p P ^ 


+ 1 


dt 


(1 - I- ) 


n 


(2.19) 


"'m 


where n , p and 6^^^ are three model parameters. According 
to Harris consistency of Bq, (2.19) with Eq. (2.17) daaonst- 
3 ?ates that Eqs. (2,14) and (2.18) are compatible with the 
model proposed by him, Ihis does not, however, prove that the 
model in Eq, (2.19) is also unique. 

In fact, it will not be an exaggeration to say that 
all the models described above result in little more than 
the curve fitting equations. Another major drawback with 
these models is the lack of any infonmation on time variation 
of the particle size distribution. In fact excluding the 

fit 

Qbarles’ model, which is^any case highly arbitra.ry, none 
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of these models give an explicit expression from which the 
particle size distributions at various grinding times may 
be predicted. It should be pointed out that a single 
parameter cannot uniquely represent the pairticle size 
distribution unless the latter exhibits a similarity or 
self-preserving character, that is the trajectory of particle 
size spectra in a batch milling operation is described by a 
single expression for the frequency function involving only 
one time dependent parameter , Ka-^ur fl3[ has presented a 
comprehensive derivation of all the best known single time 
variable parameter grinding models based on a similarity 
solution to the integro-diff erential equation of batch 
grinding jl4[ • fhis will be discussed in Obapter 2,5.3 * 

2,2 Models Based On Ihe Concept Of Discrete B^etitive 
Breakage Steps : 

Several mathematical models j 15-18, 20,21 j have 
been proposed in which the breakage process has been 
conceived in terms of discrete repetitive breakage •steps’’ 
or ’cycles’ . The concept of breakage st^ or cycle is 
easy to apply in the case of repeated shatter tests or 
single passage comminution devices such as the cone mills. 

It is, however, difficiilt to extend this concept to the 
size reduction operation in tumbling mills where due to 
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continuous rotation of the charge and the grinding media it 
is not possible to isolate a distinct r^etitire breakage 
st^ for the erfcire systan* Nevertheless, these models have 
indirectly contributed quite significantly to the present 
state of development of the grinding kinetics and , therefore, 
a brief review will not be much out of place. 

Brown [ 15 j assuned that for size reduction of coal, 
one breakage cycle can be characterized by the following 
distribution function j 


G(v,x) = 


1 - 

1 - exp (~v/m ) 


( 2 . 20 ) 


0<x<v/m, m > 1 


where x is the size of daughter particles, v is the size 
of parent particles, and m is a model parameter for the 
grinding sy^an . let, P (x) represent the fraction of 
material finer than size x after q cycles of breakage, 
then by mass balance : 

. (sx)+ /_ o(v,x) I. . (v) (2.21) 

where d.^ ("'T') = &q-1 I (2.21a) 

This model suffers from two major drawbacks - (i) 
it assimes that no particle can escape unbroken during a 
breakage cycle, which is not normally true for many grinding 
operations and (ii) the choice of the mathonatical form for 



the function 0 is quite arbitrary, which, moreover, in view 
of only a sin^e adjustable parameter in does not provide 
adequate fleKibility for adaptation to the changes in the 
process characteristics due to variation in the operating 
conditions and changes in the breakage characteristics of 
the cha.rge. It is also evident that the model is simply a 
black-box rq)resentation of the process and, therefore, 
no insight can be obtained into the breakage phenomena, 

Lpstein |l6j introduced a new concept - that of a 
'breaks ge event’, which consists of merely a single step 
or event of degradation in the breakage history of a parfeicl 
His model is based on two assumptions - (i) the probability 
of breakage of any particle during q-th step of the process 
is a constant p , 0 < p < independent of the size of 

q ""1 

the particle, the presence of other particles, and the 
number of steps that have occured prior to a given step, 

(ii) the distribution B (v, x) of the broken pieces obtained 
from the application of a single breakage event to a given 
particle of size v is independent of the size of the 
particle broken in the sense that the fraction by wei^t of 
material having size less than zv (O < z < 1) arising from 
the breakage of a unit mass of size v is independent of v, 
i.e, , B (v,x ) = B (x/v) , Thus the composite distributioii 
function arising frcm the primary breakage of a unit mass 
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of size V is given by 

(x/v) = Pq ® (xA) + (1 “ p^) A (xA) (2.22) 

where 

A(z) =0 , 0 < z < 1 

= 1 , 2 > 1 

And, the overall size distribution of the particulate 
material P (x) after g steps of breakage is related to 

H 

(x) by the equation 
^0 

P^(x) = ^ (xA) A Vi 

Incorporation of the breakage probability concept 
and introduction of the concept of breakage event are two 
most striking features of the I^stein^s model, Ihe distribu- 
tion function B which represents primary breakage of a 
particle, due to a sin^e act of degradation, provides more 
meaningful basis for process description as compared to 
the function 0(v,x) used by Brown, which represents the 
overall breakage process in which the repeated breakage of 
resulting progeny particles is also included. However, a 
priori assumption of size independence of breakage 
probability p and restriction on the primary breakage 
distribution function B should be criticized, Ebstein 
showed that under these two conditions particle size 
distribution I tends to be logarithmioo-normal asymptotically. 
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I^ora practical view point this result is not of much utility ^ 
since normally mat erials are not ground to this extent, 
j^art from this? calculation of the product size distribution 
by repeated application of Eq. (2.23) is quite a cunbersome' 
job, - especially when the number of breakage cycles is 
large. 

With a view to overcome the problem of evaluation of 
repeated integrals, Broadbent and Oallcott [17j proposed a 
practical size discretized matrix model as follows: 

Let the entire particulate assanbly be divided into 
n size intervals , where size interval i (i = 1, 2, .,,n) 
is bounded by two adjacent standard sieves of mesh openings 
and (6 = > 1) • Then, for sufficiently large 

n, the integrals, encountered in the last two models, may 
be regarded as the limiting cases of matrix multiplication, 

Bor example, if c.. . is the proportion of a typical 
particle between and x. before breakage which falls 

between and x^^ after one breakage cycle (i > j), 

then the breakage process can be represented as 

p = 0 £ (2,24) 

where p is the n x 1 frequency vector describing the size 
distribution of the product, f is the n x 1 size distribution 
vector for the feed, and 0 is the n x n breakage matrix 
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wMcli consists of n vectors representing the size distributions 
of the product due to breakage of each size interval, corres- 
ponding to one cycle of breakage. Since the size of particles 
cannot be increased by breakage , c. . = 0 for i < j s'-nd 

^ 9 J 

haice, C is a lower triangular matrix. 

Furthermore, since in each size interval i only 
-x- 

a certain fraction p- . may be selected for breakage, the 
inodsl oqus’tion beconi6S 

p = g f (2.25) 


where 

I = C P* + ( I - I*) ■ (2.26) 

'X' 

In Eq. ( 2,26 ) P is a diagonal matrix consisting of 
p . . elements and I is the identity matrix . It can be 
shown that , in general, the product of Zy cycles of 
breakage will be described by 

p = f (2.27) 


As compared to n x evaluations of integrals in Eq. 
(2.23), computation of B^ by repeated matrix multiplica- 
tion is certainly more conveniaat. Ihis is the advantage 
of size discretized matrix model over the size continuous 
ones discussed earlier. 

It should, however, be pointed out that the model 
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in Bq, (2.25) should not be considered as the discrete 
analog of %)stein*s model in Eq. (2.23). Equation (2.23) 
represents breakage product corresponding to a single event 
of primary breakage. On the other hand, the model in 
Bq. (2.25) is for one breakage cycle , in which there may 
be more than one primary breakage events i. e. , the product 
of breakage due to one event may be subjected to repeated 
breakage during the course of a breakage cycle, 

Broadbent and Callcott arbitrarily chose the 
following form for distribution function 0 i 


0 (v,x) 


■ ^1-exp (-x/vH 
^l-expC-l) " 


, 0<x<v 


( 2 . 28 ) 


and following equation was used for calculating the 
elements of the 0 matrix 

e2p(-x^/x^6°*^) - exp (-x^_^/x^6^*^) 

c. . = 

1-exp (-1) 

(2.29) 


where it has been assumed that the breakage of particles 
in size interval j can be approximated by the function 
G(v,x) corresponding to v= 3^.6 ’ , i. e. , an intemediate 

size inbetween x. and • It was further assumed that 

the fraction sheeted for breakage is independent of the 
particle size i.e, , p.- • = 7t a constant, for all i. Ihe 

Xf X 

expression for D matrix , therefore, is i 
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D=7tG + ) 1 (2*30) 

The value of % , the only unknown parameter in the model, 
was detehnined uding the standard least squares method for 
fitting datai Though, fof some grinding systems a satisfactory 
fit to the experimental data could be obtained, the v'lue of 
% showed significant variations from cycle to cycle. While 
in most of the cases the value of % constantly decreased for 
successive cycles, in some cases random variations were 
also observed. 

In order to explain the failure of the above model 
in many cases to provide an acceptable fit, - even for a 
single breakage cycle - , a new concept of 'severity of 
breakage ' was advanced by these authors |18| , They 
proposed that in certain cases the material selected for 
breakage may break more severely than that described by 
the assumed fom for the function 0(v, x), and this can be 
taken into account by raising the matrix G to a variable 
power 'ei The modified expression for D matrix can be 
written as 

1= Tt G®+(l-7t)I (2.31) 

A graphical curve-matching procedure was used for 
determination of the parameters e and n, which cannot be 
said to be very satisfactory. 
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Several comments can be made on the structure of 
the Broadbent and Gallcott^s model - (i) the choice of a 
pnrtieular form for the matrix function G is quite arbitrary. 
It is highly unlikely that a wide variety of grinding 
systans under widely different operating conditions will 
confonm to the same type of breakage pattern - differing 
only in the fraction selected for breakage . Even the 
modified form G® can generate only a fixed set of breakage 
patterns, (ii) it has been frequently observed that the 
probability of breakage varies considerably with the particle 
size • Hence, the assumption of size independence of P 
matrix is also not valid in general (iii) the model cannot 
account for the variftion of its parameters from cycle to 
cycle. It may be mentioned that Gupta and Ka-^ur |19| 
successfully simulated some of the grinding systems studied 
by Broadbent and Gallcott [ 17 j using a fixed set of parameters 
for all breakage cycles , Thus even though the grinding 
system is invariant, Broadbent and C^llcott's model gives 
a false notion that perhaps the grinding behaviour of 
particles is varying from cycle to cycle. In fact, it so 
happens that due to incorrect assumptions in the model, the 
least square estimates of the parameters ere biased to the 
particular set of grinding data used and can serve satisf- 
actorily only under the condition of nearly similar size 
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distributions. It is also obvious that the model is simply a 
data fitting scheme and is incapable of providing an unique 
identification of the grinding systan. Ihe utility of such 
models is, therefore , quite limited. 

Broadbent and CJallcott j 17 1 , Oallcott |20| and 
Lynch j 21 1 have applied the above model to ball mills also . 

It has been proposed that a fixed period of grinding can be 
considered to be equivalent to a cycle of breakage, in 
which case the model is discrete in the time va.riable . 

On the other hand the fact ranains that batch ball milling 
process is inherently continuous in time and therefore , a 
time discretized model cannot provide an adequate represent- 
ation of the process . Ihe need for time continuous model 
is felt more acutely in case of continuous mode 'of operation 
because the particles exhibit a distribution in their time 
of residence in the mill. 

It should also be pointed out that the concept of 
probability of breakage is applicable only for an infinitesimal 
time period corresponding to the occurence of a single 
primary breakage event. If the grinding time corresponding 
to a single breakage cycle is large compared to the frequency 
of occurence of primary breakage events , (for example 
2.5 min. in the work of Broadbent and Cfellcott ) an the 
particles will not break in a uniform manner. Some particles 
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will hr.ve a greater cliance for rebreakage and some less. In 
such, a case , it is evident that a single function G(v,x) 
cannot be used to describe the size distribution of the 
daughter particles . Ihis is a very serious drav/back of 
the model. 

2,5 lime and Size Continuous Particle Population 
Palance Models 
2,3.1 Mathanatical Models 

Bass, in his outstanding paper |22j , presented the 
first derivation of a phenomenological time and size continuous 
population balance model of the grinding operation for a 
description of the time variation of the size distribution 
of the particulate assembly , The model is based on statistical 
considerations and assumes that over the entire size range 
of interest, population density of particles is sufficiently 
large such that meaningful satistical averages of various 
properties associated with the grinding system may be 
defined . 

Por model representation , two functions J!(x,t) 
and G(v, x) were used, where M(x,t) dx represents the mass 
fraction of the charge in the size range (x, x+dx) and 
G(v,x) dx represents the fractional rate of formation of 

* In the original work actual mass was used instead of 
mass fraction. However, the two r^resentntions are 
essentially eg.uivalent and the basic model ronains 
unaltered. 
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particles in tlie size range (x,x+dx) from primary breakage 
of the particles in the size range (v, v+dv) . Ihree 
assumptions were made regarding the nature of the size 
reduction process - (i) smaller grains 'do not associate to 
form larger ones i. e, , there is pure disintegration of the 
particles and no agglome^^ation takes place. Hence , G(t,x) 
is defined only in the range 0< x < v and is equal to zero 
for all X > V , (ii) process of size reduction is ind<^endent 
of the constitution of the particulate assanbly and , hence, 
G(v,x) is invariant in grinding time , and (iii) the function 
G(v, x) is continuous in both the variables in the range 
0 < X < V < Xq and so is the mass density function M(x,t) , 
Detailed arguments were given in support of this assumption, 
which are quite convincing and have been widely accepted 
123 , 26 | . 

The derivation of the model equation proceeds as 
follows . Let us consider a part of the popiiLation in the 
size range (x, x+dx) equal to the mass fraction M(x,t) dx , 
it any instant of time t , the fraction leaving this size 
range is equa.l to the rate at which material finer than size 
X is being formed due to the breakage of particles in this 
infinitesimal size range , i. e. 

X 

MCz,t) dx / G(x,w) dw 
0 
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and 5 tile fraction altering this size range due to breakage 
of coarser particles is 


X, 


'0 

S fG(v,x) dx J M(v,t) dv 

X 

Hence, by mass balance it follows that 


8[M(x,t) dx] X 

= _ M(x5t) dx / G(x,w) dw 

^0 _ ° 

+ / |G(v,x} dx] M(v,t) dv (2.32a) 

X 

which simplifies to 


= _ M(xst) / G(x,w) dw 


^0 

+ / G-(v, x) M(v,t) dv 

X 


(2.32) 


This is the fundamental partial int egro-diff erential 
equation of batch grinding , 

i.n alternate independent derivation of this equation 
was lat-^r given by lilippov j23| , which was based on O- 
sophisticated probability theory approach , By similar mass 
balance considerations , ELlippov obtained 


0 - 

= f g(y) B (v,x) d^ R (v,t) 


(2.33) 


X 


and 


"^0 

= / S(t) B (y,i) 

IHl 


(2.34) 
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where R(x,t) is the statistical expected value of the 
fraction of particles coarser than size x, S(x) dt is the 
probability of breakage of particles of size x in an infinite- 
simal time interval (t, t+dt), B (v,x) is the mothanetical 
expectation of the normalized distribution function, vi/hich 
describes the size distribution of the daughter particles 
arising from one splitting of the particles of size v, and 
BCxjt) is the statistical expected value of the fraction of 
particles finer than size x. Some of the important assumptions 
made by Rilippov are - (i) breakf^ge of each particle is 
independent of its past ani the presence of other particles 
(i. e. , no particle-particle int exactions), (ii) the function 
S(x) depends only on particle size x and is bounded for the 
size interval under consideration (iii) the function B(v,x) 
depends only on v and x , and the sum of the mass of the 
daughter partiaLes , obtained from one particle throu^ 
one splitting is equal to the mass of the particle (i.e. , 

B (v,v) = 1 ) . 


It can be seen that the Bass and Rilippov models 
are time- invariant in nature • furthermore, it can readily 
be shown that both are identical in respect of the syston 
representation, too . By definition of the functions M(x,t)5 
R(x,t) and f(x,t), it follows that 


M(x,t} 


afCx,t) _ _ 3 R(x,t) 

ax ” ax 


(2.35a ) 
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aijsi 

M(x,t) dx = P(x,t) = - d^ R(x,t) (2.55b) 

Differentiating Eq. (2,34) with respect to z and making use 
of Eqs. (2,35a) and (2,35b) we obtain 

« s(x) M(x,t) - + / °s(v) M(v,t)dY 

"" (2.36) 

It clearly shows that Eqs, (2.32) and (2,36) are exactly 
similar in form. Eurther, by comparison of the Eqs. (2,32) 


(2,36) , we obtain 


S(x) = / 

0 

G(x,w) dw 

(2.37) 

S(v) 

= &(v,x) 

(2.38) 

alt ematively 



8B(v,x) 

= Glv^x) _ 

(2.39) 

0X 

V 

/ G(v,w) dw 

0 


Equations (2.37) and (2.39) also establish that both 
functions S(x) and B(t,x) in Eilippov’s model can be 
explicitly expressed in terms of Gr(v,x) , the single function 
used by Bass , Both the models are , hence , proved to be 
identical , 

Interestingly , Gaudin and Meloy j24l and Gardner 
and Austin )25| also derived the above model independently . 
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Ga3?dner and Austin derived Eq. (2,34) using functions S(x) 
and B(v,x) and named these functions as ’selection function’ 
and ’hrealcrge function’, respectively, This nomenclature 
is also in accordance with the current practice , Gaudin and 
Meloy started with number balance on particles end a 
specialized functional form for 9B (v,x)/0xj27j and obtained 
Eq, (2,36) by transfoimetion to the mass basis, Bor convenience 
in writing , the density function form of the breakage 
function 9B (v,x)/9x hereafter shall be denoted by B(v,x) . 

It should also be mentioned that application of the 
above model to tumbling mills requires an additional assumption 
, namely , the breakage process is independent of the position 
of the particle in the mill |28} , i.e. , the characteristic 
breakage event is same everywhere in the mill. However, 
if particles are fully mixed, such an assumption is not 
necessary , because in this case the average breakf ge event 
for particles of a given size ronains invariant in time , 

lastly, the above model does not take into account 
grinding by abrasion or attrition , 

2,3,2 Estimation of Model Bunctions and Model Verification 
in Continuous Size and Time 

So far , there has not been even a single report in 
the literature where the aforaaentioned model has been verified 
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k 


against the actual grinding data obtained from a physical system 


There have been two main difficiiLties in this respect , Tlrst, 
it has not been possible to determine the selection and 
breakage functions S(x) and B(v,x) from purely theoretical 
considerations based only on the fundamental physical 
properties and dynamics of the mill-material system . Second, 
there are several physical limitations which preclude the 
direct detemination of these functions using experimental 
methods. It is virtually impossible in practice to isolate 
the material of size x exactly . Equally difficult is the 
task of performing a grinding experiment with tumbling mills 
in v;hich no particle undergoes more than one i^litting , 

] hese difficulties do not permit us to take advantage of the 
well known methods of syston analysis such as the tre.cer 
techniques . A brief account and criticism of the work done in 
this area is given below . 

Bass \ 22 \ subdivided the size range 0<x< x^ into n 
size intervals of very narrow width such that it is reasonable 
to assume that the fraction of the product of breakage 
ronaining in the original size interval is negligible , Under 


this condition he proposed a design of experiments which will 


permit determination of various parameters v^^ , where 

\-i ^1-1 


^-1 - ^ 


/ G(v,x)dx dv 



( 2 . 40 ) 



I 
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and and are upper and lower size limits, respectively, 

for i-th size interval. It was further proposed that a 
generalized polynomial in v and x could he used to approximate 
the function G(v, x) and the coefficients of this polynomial 
could he determined hy minimizing the following error 
function via the method of least squares 


n-1 n 

er = 2 Z C ■^ik 1 - 

k=1 i=k+1 IK K-T K 


\-1 ^ i -1 2 

/ / G(v,x)dx dv 3 

\ ^i 


= + w* G(x,x) (2.41) 

where er is the error function and w is a weightage factor. 

In principle the method appears to he sound. However, 
there are several draw-hacks which make it highly inpracticahl e 
- (i) it requires that extremely narrow size intervals 
should he chosen , As to^shown later in Chapter 6 , size 
intervals as narrow as corresponding to the ratio 6=x. ./x. 

I il. 

s= 4^2 paoinot satisfy the condition of ne^igihle inner- 

hreakage . Therefore a large numher]^^x(n-1)/2 v’s must he 

determined hy performing a equally large number of experiments 

, (ii) due to the presence of unavoidable experimental errors 

in the measuranent of the size distributions, v. . values 

J 

will also be in error. In absence of any knowledge of the 
magnitude of errors introduced in v's due to the experimental 
end conputational procedures, the least squares method 
does not guarantee a real identification of the system, and 
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(iii) as to be shown later in Chapter 8 , grinding by abrasion 
cannot be ne^ected and , hence , Bass's assumption of Gr(x,x)=0 
is not valid , in general . Ehe fomulation of the error 
function in Eq. (2.41) is also, therefore , not quite valid , 
lastly , no daaonstration showing the potentialities of the 
above identification scheme using any real grinding data has 
been made . 

Based on some theoretical consideorations Gaudin and 
Meloy j24 , 27 j derived the functional foms for both the 
functions S(x) and B(v,x) , Ihey did not, however , describe 
how the various contents in these expressions were to be 
determined for a given grinding system < Assigning arbitrary 
values to these parameters they computed size distribution 
curves for various grinding times and concluded that the 
nature of these curves was very much similar to those obtained 
in the actual practice . Again , the model was not tested . 
against the grinding data any real system . 

Gardner and Austin [25j proposed a procedure whereby 
the size discretized analogs of selection and breakage 
functions could be used to determine the functions S(x) and 

* It may be pointed out that this condition means that 
B(x,x)=0 for all x , which in view of various findings does 
not seen to be true . 
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B(v,x) . These authors assumed that the selection function 
S(x) for X = (x^ + Xj^_^)/2 is same as the selection parameter 
fp.r i-th discrete size interval , and the cumulative 
breakage function B(v,x^ ) for v = (x^ + x:^ same as 

B. . , It may be mentioned here only briefly that both S. 

X j J 1 

and B. . are determined by S(x) and B(v,x) jointly , and , 
m > J 

therefore , exo^t in special cases , as to be discussed 

later in Chapters 5 and 6, it is impossible to find out S(x) 

from S. or B(v,x) from B. , alone , Hence , this procedure is 
1 1)3 

completely arbitrary arisi::a:hsiD 2 d , and has no relevs.noe 
whatsoever to the determination of the functions S(x) and 
B(v,x) , 

In another paper j 29 ) Gardner and Austin used smooth 
curves of S(x) and B(v,x) to determine the size discretized 
parameters using the following two equations i 

= [ S(x^) + S(x^_^) ]/2 (2.42) 

\,j ^ + i (x^ , x^__^)]/2 (2.42a) 

It should be pointed out that in addition to the 
fact that Eqs. (2.42) and (2,42a)are totally incorrect and- 
, it hardly makes any sense to convert discrete 


* The definitions of these parameters are given in 


Chapter 2,4 
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parameters to continuous functions and then back convert 
continuous functions to discrete parameters . Moreover , it 
can be proved that the process of conversion used by them is 
not reversible . 

In conclusion then it can be said that the population 
balance grinding model in continuous size and time has not 
been properly ^rerified . However , in the next section we 
shall witness some indirect evidences in suppoit of the 
model 5 ?/hich have been obtained using intuitively appealing 
functional foms for the functions S(x) and B(v,x) and 
comparing the experimental grinding data with the exact 
solution to the basic integiXD differential batch grinding 
equation , 

2. 3. 3 Some Analytical Solutions to the Population Balance 
Batch Grrinding Equation in Cfontinuous Size and lime 
The only general solution to Eq. (2.36) , the 
fundamental int egro-diff erential equation of batch grinding , 
is due to Bass |22l , who obtained the following solution in 
terns of iterated kernels s 

X 

M(x,t) = f(x,t) exp f-t / &(x,w) dw 1 (2.43) 

0 

where 

^0 t °° 

f(x,t) = M(x,0)+ / dv / S H. (x,tjv,t') M(v,0) dt' 

X 0 i=1 ^ 

(2.44) 





(2.45) 


H 


1 


IL (x,-fc| v,t') = 

and 

G'(v,x) = G(v,x) 


/ dw / G(w,x) H. y,t*) dt 

V t f 

( 2 . 46 } 

V X 

eap [ - t'( / G(v,u) du - / G(x,u) du)] 
0 0 

(2., 47) 


Because this solution is quite ccmplicated and computationally 
quite time consuming , it has not found any practical use' 
thus far , 

Some simple solutions have been derived for specialized 
functional forms of the functions S(x) ard. B(v, x) . Gaudin 
and Meloy |24| showed that for 


S(x) = A X 

(2.48) 

and 


B(t,x) = x/v 

(2.49) 

the solution to Eq, (2,36) in cumulative retained 

mode is 

R.(x,t) = R(x,0) exp (-A x t ) 

(2.50) 

where A is a constant , 


Austin et al, |30j presented a more general solution 

assuming the following restrictive relationship 


S(v) B(v,x) - A g(x) 

(2.51) 

Since , B(x,x) = 1 , it follows that 


S(x) = A g (x) 

(2.52) 



and 


i(v,.) = ^ 


(2.53) 


j I 


Jurtheimore , since B(t,x) is a monotonically non- decreasing 
function .,in x , g(x) should also satisfy this condition.. 

And , since B(v, 0) = 0 , therefore , g(0) = 0, It follows 
that S(0) = 0 and S(x) is also a monotonically non-decreasing 
function of x , The solution to the grinding equation is 

R(xjt) = R(x,0) exp {J - A g(x) t ] (2.54) 

0^ Xq 

M(x,t) = esp [ -A g(x) t ] |“M(x, 0) + At / M(v,t)dvj 

(2.55) 


When g(x) = x, Eq. (2.54) reduces to the solution of Gaudin 
and Meloy j24| , g(x)= - In g(x) gives the solution of 
Iferris j3lj , and , for g(x) = x®, it gives the well known 
Rosinr-ESmmler size distribution equation . It should be 
mentioned that for many grinding systems studied in the 
past jl7, 13t 32-36 j , Eq. (2,54) was found to hold 
reasonably satisfactorily , This may be taken as one indirect 
evidence in support of the grinding model. 

Using the functional forms 
S(x) = A x“ (2.56) 

B(v,x) = ( I (2.57) 


King j26| obtained the following analytical solution in 

A. 




, \ 
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cumulative fraction finer mode 

J’(x,t) = ( hC-p/a , p/a ; 1; - A x “t , - A x“t) 

0 ^ 

^0 P 

+ / [( ~ 5'(w,0) - M(w,0)) ( f exp (-Aw“t) 

(p/a , 1, (w“ - x°^) A t) j dw (2,58) 

where h (a, b, 1, x, v ) is the generalized hyper geometric 

series in two variables and (a, 1, z) is the confluent 

hypergeometric function , Phou^ , this restricted solution 

to the grinding equation is computationally more convenient 

as compared to the genei®.! solution obtained by Bass , however, 

generation of the values for the hypergeometric series and 

confluent hypergeometric function on digital computers is 

quite a problem , Por hand calculations reference can be made 

to the tables [37 1 and the graphs given by Zing j26| , King 

gave an analytical solution for piecewise continuous selection 

function also, She selection function exhibits a maxima at the 

* 

point of discontinuity x and is proportional to x for 
X < X ' (a^ > 0) and to x for x > x (a^ > O) • 

Under the restrictions assumed by King and given in 
Eqs. (2.56) and (2,57) , Kapur 1 14] obtained a very simple 
and useful similarity solution to the grinding equation 
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where 



a 

pTo 

(ati) r(p/a) 


_ 1 1 “ }' (^ /a) 

“ rC(^+l)/a ] 


( 2 . 60 ) 


(2o6l) 


and r is the gamiiia function . The scale factor , 

the first moment of the distribution M(x,t) , varies with 
time in the following manner 


(t) 

dt 


a+1 

Ah [ ] 


( 2 , 62 ) 


Kapur jl3-14j bas shown that ylL^(t) is proportional to 80 per 
-cent passing size , It can be shown that it is also proportional 
to X as defined in Bq. (2.4-) fAppendix 1] <. Moreover , for 
a grinding period of moderate duration y it is frequently a 
fair approximation to assume that the power draught is 
constant |38j , and therefore , E , the g 2 ?oss energy input, 

s 

is directly proportional to the grinding time t « Thus 
replacing by x , (a+1) by r ard t by , Eq. (2.62) 
can now be written as 

i- 

dE 

— £ = _ (2.62a) 

dx X 

where Cg is a proportionality constant , v/hich incorporates 
the constants A and h of Eq. (2,62) « It will be seen that 
Eq. (2,62a) is similar to Eq. (2,1) , with the difference that 
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instead of E now E^ is the energy variable . As pointed oxrb 
earlier in Chapter 2.1 , for a verifi-'tion of the energy - 
size reduction laws it was invariably assumed that E or E' 
are directly proportional to E , an assumption for which 
no rational basis exists . On the other hand , Eq. (2.62a) , 
derived from considerations of particle population balance 
has in fact , exactly the same form in which all energy - size 
reduction laws have actually been verified . It is interesting 
to note that according to Kapur’s derivation , the constant r 
is in no way related to any particular energy-size redtictTion 
law , Instead , it is function of the exponent a , which 
refers to altogether different phenomena - those which 
detemine the size dependence of the probability of selection 
for breakage of particles . It should , hence , be wry clear 
that in Eq. (2.62a) r = 2 corresponding to a = 1 does not 
necessarily imply that Rittinger's law (i.e. , AE^aq) is 
valid (and similarly for r=1 and r=1.5) . On the contrary 5 
two mill-material grinding systems exhibiting same value of 
a - the possibility of one obeying Kick’s energy law and the 
other Rittinger’s law cannot be ruled out , 

Two important points emerge from above discussion . 
Eirst , the extensive literature which is devoted to the 
verification of Eq. (2.62a) in one form or another , gives 
a wide support to the population balance model used by Kapur . 
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mill j39j were shown to be approximately self-preserving over 
a fairly wide range of the values of />i^^(t) , This may be taken 
as another indirect evidence in support of the grinding model .* 
Kapur [14| also indicated that the crude estimates of the 
parameters a and p could be obtained by simply matching the 
experimental similarity distributions with the standared 
theoretical ones for various values of a and p , Austin et al. 
|50j also suggested a graphical curve matching procedure for 
determination of these two parameters . Interestingly, neither 
Kapur nor Austin et al, estimated these parameters for any 
real grinding system , Therefore , the direct verification of 
the model has not been registered , thus far . 

It may also be mentioned that under the constraints 
of Eqs, (2,56) and (2.57), exact solutions to the grinding 
equation have been derived by Tilippov [23| and King |26[ for 
the initial condition of an inipulse situated at size Xq . This 
result is of theoretical in 5 )orfeance only , because in practice 
it is not possible to experiment with a material of exactly 
a single size . 

2,3.4 Time Modified Pirst Order Kinetics 

iferris {38| and Kapur [13,14| have presented 
solutions to the fundamental integro-diff erential equation of 
batch grinding assuming specialized time dependent functional 
forms for the selection function . Harris assumed the following 
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restriotiye forms for tlie two basic functions 

• S(x) = - At^ In g(x) (2.63) 

and 

B(v,x} = (2.64) 

In g(v) 

where the constant ^a’ serves the piurpose of providing small 
systonatic changes in the selection function v/ith the grinding 
time t . (Ehe solution to the grinding equation Eq. (2, ‘34) is 

P(x,t) =1-1 1-E(x,0) ] [g (x) ] (2,65) 

Iferris noticed that when a 0 j Eq. (2.65) does not give 
consistent results . Por example , if the grinding time 
interval t is divided into two subintervals of durations t^ 
and ±2 » respectively (such that t = t^ + t 2 ),and P(x, t) is 
formed by taking the product from first interval as the feed 
for the second interval , the result obtained is not the 
same as that obtained by directly substituting t = t'^ + tg 
in the Eq. (2,65) , Harris , therefore , inferred that the 
constant 'A’ should in fact be a function of E(x, 0) such 
that consistent results are obtained irrespective of the 
method of calculation used . He did not , however , make any 
attaapt to find out this relationship , 
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Eapur 1 13, 14 [ assumed the following form for the 
selection function 

S (x) = A x“ 1 (t) (2.66) 

where 1 (t) is an unspecified function in grinding time t , 
which represents the influ^xce of particles' environment on 
the selection process . If the breakage function is governed 
by Eq. (2,57) , then the first moment of the particle size 
distribution varies with grinding time as follows 

TTtF^t = - -k h [/-^^(t) ] (2.67) 

and a similarity solution to the int egro-dif ferential 
equation of batch grinding is given by Eqs. (2.59) , (2.60), 
(2.61) and (2,67) . In particular , for derivation of some 
energy - size reduction relationships Kapur |13,40| assumed 
the following form for function I 

I (t) = 1/ A^(t)^ (2.68) 

where a is a constant 

It should be pointed out that the grinding time 
itself does not directly influence the grinding operation 
in batch tumbling mills and , therefore , the grinding system 
as a whole cannot exhibit any explicit time dependence . 
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Any apparent time variability in the behaviour of the grinding 
system may be due to several reasons . One of the possible 
reasons , acco 2 ?ding to Austin et al. [41 [ , is due to the 
interaction between the distributions in strength of particles 
and in force applied by the grinding media . I-tilca et al. |42| 
and Kapur (13,40| have indicated the possibility of the 
influence of particle size distribution on the rate of breakage 
of particles . Also, Harris |38| has discussed how the 
proportion of active grinding sites may change with increasing 
fineness of the material . It will thus, be seen that the 
time variation of the breakage characteristics of the particles 
is determined by the initial conditions such as the strength 
distribution of particles at t = 0 , the particle size 
distribution at t = 0 , and several other similar possible 
factors [281 . Since , in different experiments, these 
initial conditions need not be the same , it is obvious that 
the time variability cannot be described as a unique function 
of time for all grinding experiments . In this respect , 

Kapur's approach should be less restricted since he used the 
first moment of particle size distribution to describe the 
time variation of the selection function . The 'a priori' 
assimnption of the particular functional form in Eq. (2.68) 
may , however , be questioned . 
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Ifekajima and Tanaka { 43 1 derived the following 
approximate solution to the batch giinding equation under 
the restrictions given in Eqs. ( 2 . 56 ) and ( 2 . 57 ) ! 

E (x,t) 1- [ 1 - E (x,0) ] exp f- (a^ A x“ t)^^ j ^ 2 . 69 ) 

where a^ ard a^ are constants determinable from p/a as 
described in their paper . These authors have tried to show 
that Bq. (2,69) , which is virtually the same as Harris' 
equation Eq. (2.65) for the following form of function g (x)|38| 

®3 ^4 

e M = [ 1 - (- 5 - ) 1 ( 2 . 70 ) 

can be derived without assuming a time dependent selection 
function . In a recent article |44| Austin et al, have pointed 
out that the derivation of Eq. (2.69) is suspect because the 
constant A is not a function of E(x, 0) and , hence , it 
leads to the same paradox which was noted by Harris 1 38 | 
in the case of Eq. (2.65) . 

Before concluding this section a few comments are 
necessary regarding the need and justification for the use 
of time modified first order disappearance kinetics . It will 
be seen that if the selection and breakage functions are 
interrelated through a common function in the manner given 
in Eqs. (2,52) and (2.53) ? then according to Eq. (2.54) a 
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In In |_ R(xjO)/R(x,t ) ] / In t plot (jlLyavdin - WeibiiLl plot 
|45-47l) should result in a straight line of slope unity , 
There is a general opinion that as selection function and 
breakage function refer to altogether different phenomena , 
such an interrelationship between them is highly suspect , or 
even unacceptable 130,511 . It, therefore , follows that the 
deviations from a slope of unity as observed by ^rris 
j45-46j are orOLy natural and to be expected . But , at the 
same time , unless we are able to filter out the effects due 
to the violation of the above mentioned interrelationship 
between the two functions , there is hardly any need to 
assume a time d^endent selection function in the first place. 
In fact, at present , it does not seem possible to aralyze 
any inherent time dependent effects properly mainly because 
we do not have an adequate basis to obtain a precise knowledge 
of the selection and breakage functions * 

In conclusion , the time modified first order 
kinetics as introduced by Harris through Bq. (2,64) has no 
physical basis wfee^«©=«6©^ , Moreover , the graphical 
demonstrations of a wide applicability of Bq. (2.65) by 
Harris |45-46| should also be accepted with some reseivations 
as there is an inherent tendency for linearization in the 
plotting procedure based on highly compressed log log versus 
log scales « 
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2.4 Size Discretized Time Oontinuous Particle Population 
Dalance Models 

As it can clearly be seen from tbe discussion in 
tile previous section , because of the problems in the 
identification of the grinding system , as v/ell as , the 
complexity of the general solution to the fundamental 
int egro-diff erential equation of batch grinding , the size 
continuous model has found little use in practice . On the 
other hand , since the particle size distribution is invariably 
determined by sieve analysis and is available only in frequency 
form as mass fractions in discrete size intervals , there has 
existed a valid motivation for an appropriate practical size 
discretized model of batch grinding kinetics , first proposed 
by Sedlatschek and Bass |48( in the year 1953 . Since then , 
many authors 1 22, 28, 42, 48 - 54 1 have discussed the 
rational for and justification of the discrete size mathematic- 
al models formulated by them . A critical review of these 
models is presented below . 

The notation used in this discussion is as follows . 

A discrete size interval is defined as the size range 
bounded by the size of the mesh openings of two adjacent ^ 
sieves in a standard series , normally characterized by a 
constant sieve size ratio 6 (> 1) . The entire size range of 
interest is divided into n size intazrvals . Por i-th size 
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interval tiie size of "tiie mesii openings of the upper sieve is 
denoted by ^ and that of the lov/er one by x^ . The upper - 
most size interval is represented by i = 1 and the maximum 
size limit of the particles by x = Xq , The mass fraction of 
particulate solids in i-th size interval at grinding time t 
is denoted by (t) . In the literature wherever the actual 
mass of the particles has been used in the model formulation 
it has been replaced , for convenience and uniformity , by 
the^fraction in the following discussion . This is permissible , 
because the basic treatment remains unaffected . 


2,4.1 Mat hemat i cal Mo del s 

i^sed on statistical considerations Sediatschek and 
Bass |48| proposed that the disappearance kinetics of particles 
in a small size interval i can be described by the following 


equat ion 

d M^ (t ) 



M^ (t) 


(2,71) 


where v. . is the specific velocity of breakage of the 
JL ^ 3. 

particles in the i-th size interval . This was confirmed by 
plotting In M^ (t) vs . t , which resulted in straight lines 
with only slight deviations . It may be noted that the slope 
of this straight line is equal to the numerical value of the 
constant v^ ^ . further , a specific rate of appearance v^^^ 
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wss defined as the fractional rate of formation of material 
in the i-th size interval due to the breakage of the material 
in the g-th size interval (j < i) . I^om mass balance it 
follows that 



( 2 . 72 ) 


In general , the overall mass balance for any size inter'/al i 
can be written as 


d M. (t) 
dt 


i-1 

= - V. . M. (t) + Z 

1,1 1 


^i,D 


(t) 


, i = 2, 3, ... (2.75) 

Urom the simultaneous solution to the set of equations given 
by Eq. (2,73) and a series of grinding experiments carried 
out by them , Sedlatschek and Bass obtained the values of 
the coefficients v. . 's for alumina - ball mill grinding system, 
which was characterized by four size intervals . These values 
produced a satisfactory simulation of the actual grinding 
experiments . 

In a later publication Bass j22| presented the 
derivation, of the size discretized model in Eq. (2.73) 
starting from the fundamental int egro-diff erential equation 
of batch grinding ( Eq. 2.32 ) . The derivation \v8-s based 
on the following assu 2 i 5 )tion 

^1-1 - << t^i-1 + 


(2.74) 
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such, that the following two conditions are satisfied - (i) 
the mass of particles remaining in the original size interval 
after undergoing breakage is negligible ; and (ii) in detemin- 
ation of the following average values 


""i-1 

/ M (xst) dz = M(x^,t)[ - x^] (2.75) 

Xi 

^ i -1 = ^ 0-1 " i -1 

/ dx / Gr(v,x) M(v,t) dv = M (x. ,t) / dx / G(v 5 x)dv 

X. X. ^ X. X. 

0 1 0 1 

(2.7f) 

it is valid to write 

M ( , t) = M , t ) (2.77) 


The exact mathonatical definition of the parameters v. . 

i»0 

is then given by Eq. (2.40) . In Ghapter 2.5 the rel event 
criticism of this model in connexion with the first condition 
of Bass , has already been given , which need not be repeated 
here . A more practical and mathenatically less stringent 
condition for time invariance of model parameters shall be 
discussed in the main body of this wo]ck . 

Reid 1 49 1 proposed a 'practical batch grinding 
equation ' based on ’practical parameters' , which can be 
directly detemined by experiments* . He noted that those 


l.l.T. K/'^FUR 


fBRARY 


No* 


461S9. 
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particles , which after undergoing breakage remain in the 

original size interval , cannot be measured experimentally , 

Therefore , he considered it necessary to redefine the breakage 

is 

event such that a particle is broken only if it^no longer 
retained in the original size interval , Reid pointed out 
that due to the new definition of breakage , the breakage 
function will vary with the screen ratio 6 , and defined a 
modified breakage function B (w, x, 6) , which gives the 
size distribution of the products of primary breakage of 
material originally in the size range w to 6w . Also, a set 
of breakage parameters B. . were defined as - the weight 

J 

fraction of the naterial broken finer than size x. which 

D 

reports to size interval i , when the material of size interval 
3 undergoes prinBry breakage . Another set of parameters , 
fractional ra.te of breakage for size interval i , v/as defined 
as 

= - dM^ (t)/M^ (t) dt (2.78) 

Reid pointed out that also incorporates the above 
^definition of breakage . Purther , he assumed it to be 
Independent of time , 

. The particle population balance , when expressed 
in terms of the mass fraction and the parameters and 
B. . , leads to the following kinetic expression for the 
batch grinding operation 
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dM^(t) i.l 

dt % M^(t) , i = 2, 3, ... n 

(2.79) 

Reid mentioned that the time independent breakage 
parameter B reflects the inherent assumption that the size 
distribution within a size interval remains constant , If 
sieve size ratio 6 is small , natural size distributions can 
be well approximated by straight lines over single size 
intervals the assun^tion of a constant size distribution 
is valid , He , however , failed to realize that the same 
argument can be used to justify the use of time- invariant rate 
parameter in Eq. (2,78) . Instead , he first went in a 
roundabout way to show that in general is not independent 
of time ( interestingly , even here his basic approach to 
the ai-guments is not correct ) , Then, he invoked some 
experimental results on the decay of single size coal |25|, 
alumina j48j and quartz [ 50 | ground in batch ball mills and 
pointed out that the time invariant rate parameter in Eq, (2.78) 
could provide a close approximation to the actual decay kinetics 
over the time range investigated , and hence, the assumption 
is valid , 

The major criticism of Reid's work is that he 
formulated the size discretized model in isolation from the 
more fundamental size continuous model , and failed to provide 
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a meaningful interpretation to S. and B. . parameters in 

terms of tiie basic functions S(x) and S (v^x) , The lack of 

insight into the structure of these paziameters prevented him ' 

from obtaining the exact mathematioai conditions for time 

invariance of the two sets of parameters . As to be shown 

later in this work , these pajcameters are uniquely inters 

related , but Reid's work gives a false inpression that the 

two sets of pai*ameters can assume independent values . 

RLnally , it is very clear that both the definition of the 

parameter as given in Eq, (2.78) , as well as, the definition 

of the breakage parameters B. . given by Reid , do not refer 

J 

to the particles remaining in the original size interval , It 
is then quite unnecessary to modify the definition of the 
breakage event . Moreover , it can be shown that in certain 
cases it will lead to confusing and misleading cont2?aidi ctions, 
Ror exanple , if chipping and abrasion are the dominant modes 
of grinding for particles in a particular size range , such 
that practically all the broken particles fall back in the 
original size interval , then , according to Reid's definition 
of breakage , even though the particles do not break , 
nonetheless the material gets ground , which is absurd , It 
must be emphasized that even those broken particles which 
do not break according to Reid's definition , do contribute 
to the values of the rate and breakage parameters . 
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A con^jarison of Eq. (2,71) v/ith Eq. (2.78) and 
Eq. (2.73) with Eq. (2.79) reveals that the two models due to 
Bass and Reid ere basically indentical , and it follov;s that 

- ^i,i (2.80) 

and 


B. . 


V. . 


(2.81) 


We also define a cumulative breakage parameter 
B. . as follows : 



S 

k>i 



( 2 . 82 ) 


The mathematical model of Mika et al, |42| , which 

can also be represented by Eq. (2.79) t is based on the 

assumption that the fraction of broken particles which remain 

in the original size interval is negligible i.e. , B. . = 0 , 

1 5 1 

As to be shown later in this work (Chap. 6) , this assueption 
is grossly incorrect for practical sieve size ratios of 
and ^2 , Oonsequently , the definition of B. . , namely, 
the proportion of the primary breakage product of the j-th 
size interval v\^hich reports to the i-th size interval , is 
evidently incorrect , Eor the same reason , the use of term 
'Selection Eunction ' for parameters is also incorrect , 
because all the particles select ed&r breakage do not fall 
out of the origiiBl size interval . It should be emphasized 
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that the experimental measurements of the breakage pa2?a.m6ter 
do not r^ resent the actual size distribution of all the 
progeny particles , rather represent the distribution of 
only a fraction of such paiticles which fail out of the original 
size interval , like Reid , these authors also did not make 
any attempt to identify the size discretized parameters in 
terms of the basic functions S(x) and B(v,x) , They assumed 
that the breakage parameters -B.. . are time independent while , 

X, J 

the time variability of parameters was attributed to the 
influence of size distribution of the environment on the 
breakage probability of the pairticles . 

Preeh et al, |51 , 52 1 made an ad-hoc stipulation 
that the ground material is fed forward to only next three 
or four smaller size intervals . This was perhaps motivated 
by the resulting convenience in the system identification 
due to a reduction in the number of unknown parameters , They 
called the parameters B. . as ’inter-size flow coefficients’, 
which in light of the discussion above seems to be a more 
appropriate term than the commonly used 'breakage function' , 

The arbitrary nature of this model should be criticized 
because mere data fitting is hardly equivalent to. a realistic 
modelling , 
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linally , it may be mentioned that the formal 
mathanatical structure of the grinding models used by various 
investigators , such as Kelsall and co-workers j 55 > 56| , 

Reid and Stev\/art 1 57 1 , Herbst and co-workers |58 - 60j , 

Klimpel and Austin | 61 1 , Gardner and Sukanjanjtee |62, 63| , 
Moreira et al , | 64| , Kapur and Agrawal [ 65 - 67 1 and many 
others , is adequatd.y represented by Eq. (2.79) . iHie ; 

differences are only in the terminology used and the 
definitions assigned to the parameters S- and E. . , some of 

1 1 5 0 

which have already been discussed earlier . 

An ircportant point for discussion is the precise 
relationship between the discrete parameters S. and B. . , 

^ 1 5 J 

and the size continuous functions , S(x) and BCv, x) . With 

the exertion of a recent note published by Austin and 

Ehatia |68| , there seems to be a complete void in the 

grinding literature on this subject . In fact , a number of 

investigators , without appreciating the implications ^involved.,, 

have used these terms interchangeably [ 69 - 7 1 1 . One fails 

to understand from the literatuie that if s(x) •:<: x^ , then, 

cx 

how 3. is also proportional to x. , and , similady , if 
B(v, x) = (x/v)^ , then how B. . = (x. /x^) j 30 , 62 , 72 j . 
Examples can be quoted from the literature l62| where the 
functional forms used for discrete parameters are not consis- 
tent with the assumed functional forms for continuous 

f# 


i 
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functions and , again , where the form chosen for selection 
parameters is inconsistent with the one chosen for breakage 
parameters | 62 , 72 - 75| . In the latter case , it is 
iii 5 )Ossible to find two appropriate mathanatical functions 
for S(x) raid B(v,x) , which when used in the expressions for 
parameters S. and B. . shall give rise to the assumed express- 

-L 1 5 J 

ions for the two sets of parameters . An erroneous procedure 
of interconversion by Gardner and Austin |25 - 29 | was 
pointed out in a previous section johap. 2.5 [ . 

Assuming reasonably small size intervals , such 
that , over the size range x^ < x < x^ ^ the particle density 
M(x,t) can be approximated by M^(t)/ (x^ "I ~ ^i ^ from 

the very physical definition of the parameter ? i.e, , the 
fractional rate at which material is falling below the bottom 
size of the interval i due to the breakage in single size i , 
Austin and Bhatia j68| obtained the following expression 


X 


S. 

1 


/ 

X. 

1 


i-1 


S(v) B(v,x^) dv 


^i-1 -^i 


(2,83) 


In the same fashion they obtained 
/ sCv) B(v, Xj^) dv 

= i 


b-i 


X. 

D 


S.. B. . 
d d 


(2.84) 
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from which the expression for B. . follows immediately. 

In the present work the exact definitions of S. and B. . 

1 1,3 

parameters have been derived in Gh^ter 5 , Eqs. (2,83) and 
(2.84) are shown to be a specialized case of the exact 
r el at ion ship s . 

Some investigators j 53-54 1 have used a different form 

for representation of the discrete size grinding model. Ihe 

grinding system is characterized by two new sets of parameters- 

the selection parameters Sj^ and the distribution parameters 

b. . , which by definition are exact discrete size analogs 

of size continuous selection and breakage functions S(x) and 

B(v,x) , Ihe selection parameter s^ is defined as the 

fractional rate at which the material in size interval i 

gets selected for breakage • Ihe distribution paraneter b. . 

r^resents the fraction of all material selected for breakage 

in the size inteirval j , which directly reports to the size 

inteiTval i without rebreakage , I&ilike the breakage parameters 

B. . , which refer to only those particles which break out 
i» 3 

of the size interval 3 , the distribution parameters b. . 

^ J J 

give the true discrete size breakage distribution and 

therefore , b. . is also a valid member of the set of the 
3 >3 

distribution parameters . The equation describing the 
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the disappearance kinetics of a sin^e size fraction can now 
be written as 


dli. (t) 

■"'”t"'" = “ + Sj_(t) b^^^(t) M^(t) (2.85) 

-And, the batch grinding equation is 


mAt) i 

-" I t " " ^ "" 

, i = 1s2,..n (2.86) 

It should be noted that in contrast to Eq. (2.79) the upper 
limit of suirmation is now i instead of i - 1 .A conparison 
of Eq. (2,78) with Eq. (2.85) and Bq. (2.79) with Eq.(2.86) 
gives 

Si(t) = Sj_(t) ( 1 - \^i(t) ) (2.87) 

and 


S.(t) B .(t) = s.(t) b. .(t) 

J J 

Ecom Eqs. (2,87) and (2.88) it follows that 


( 2 . 88 ) 








(2.89) 


It will be seen that s- , b. . and b. . cannot be 

J d y J 

explicitly expressed in terms of S. and B. . alone , Hence, 

0 J 

only one way conversion given in Eqs. (2,87) and (2.89) is 
possible . It should be pointed out that from the practical 
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point of Tiew there is hardly any reason to prefer Eq, (2^86) 
oyer Eq. (2.79) . This is because s- aid b. . can always be 

X X ^ X 

lunped into one parameter S. and the product s- b. .is 
also replaceable by S. B. . , In other words, since b. . 

J X 5 J I5 X 

is not required for the formulation of the model separately 

, it is unnecessary to have n more model parameters . .. .Ari: the 

same time , it should be admitted that the parameters s- , 

D 

b. . and b. . provide more insight into the brealcage 

j , J X, j 

behaviour of the pairbicles than the rate and breakage 

parameters S- and B. . . Eaxim the point of view of theoretica]. 
J X, J 

analysis, therefore, the system model representation in 

Eq. (2,86) may prove to be more useful. 

Olsen 1 53 1 and Whiten [54j have derived analytical 

expressions for s. and b. . in terms of the functions S(x:) 

X 1,3 

and B(v,x:) , Olsen rightly obtained the following expression 


for s^ 


^i-1 


M(v,t) S(v) dv 

P 


He assumed that B(v,x) = (x/v)' and erroneously concluded 
that it implied b. . = b, . . His expressions for the 

X;rJ 1— J 

distribution parameters are 5 


( 2 . 90 ) 


h,i = 


(5)^ 

p(i-D) 




(2.91) 

(2.92) 
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Tile relationships given by Whiten are : 

A-1 

/ S(v) dv 

3. = _i 

^i-1 “ ^i 


(2.93) 




^ 0-1 ^ i-1 

/ / B(Vi,x) dx dv 

^i 

"^3-1 ■ ""d 


(2.94) 


It will be shown later in Chapter 5 that Eqs. (2,91) and (2,92) 
as well as Eqs. (2,94) and (2,93) all are incorrect , Moreover, 
we shall show that these parameters cannot be assigned unique 
values , because they are functions of the particle size 
distribution withi.-the discrete size inteirvals. 

Henceforth, we shall refer to grinding model in 
Bq. (2.79) as model 'A* and the one in Eq. (2.86) as model ’B*. 


2.4.2 Analytical Solutions to the Size Discretized Batch 
Grinding Equations 


Beid j49j was first to obtain an analytical 
solution to the set of n differential equations given by 
Eq. (2,79) • jE!he solution for un-equal rate parameters is 


as follows : 

M^(t) 


1 

2 


3-1 


h. . exp (- S.t) 

i > d ” d 


j 


(2.95) 
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where 




i-1 

S 

k=0 


S, B. 
k 1 


»k 


S. - S. 

1 3 


k i 


( 2 . 96 ) 


i-1 


\,1 = \,k ■ i 1 (2-9'?) 

and lij^^=M^(0) (2,98) 

In case two or more size intervals may exhibit same 
value of the rate parameter, i.e. , , then as shown 

by Olsen 1 53a] , Eqs, (2,96) and (2.97) should be modified as 

i-1 

(2.99) 


h, 


r,D 


B. 


L ^ " i,k --:-,3 


k=3 


and 


i-1 


h. . — M.(0) — s h. , i 7^ 1 

1,1 1 . 1,K 

k?^3 


(2.100) 


The set of n simultaneous differential equation in 
Eq. (2,79) is conveniently represented by a matrix equation 


dM(t) 

= - 1 ; “ 5 I I (t) 


( 2 . 101 ) 


where M(t) is a n x 1 vector of mass fractions M^(t) , B 
is a n X n lower triangular matrix of breakage parameters B. . 
and S is a diagonal matrix of rate parameters . 1. matrix 

solution to Eq, (2.101) due to Herbst and Mika J76j is 

-1 

M(t) « I V(t) N M(0) 


( 2 . 102 ) 
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or 

M(t) = 0(t) M(0) 


(2.103) 


where the elements of matrices N and J are , respectively 

(2.104) 






0 

I < D 

i-1 \,k \ 

i = 3 


i > d 

e:p ( - S^t) 

i = 3 

0 

i 3 


(2.105) 


Ihe solution in Eq, (2,102) in the form of an input-output 
relationship is of great importance • fhe elontnt o. . of 
the input-output matrix Q(t) directly gives us the fraction 
of matrial originally in size interval j which after breeikage 
for time period t is found in size interval i , It is interes- 
ting to coup are the matrix 0 with Callcott*s mil i matrix E , 
Needless to say that 0 being esplicitly defined in terms of 
physically identif icable quantities S. and B. . , and being 
a continuous function of time , is more meaningful and useful 


than E • 

It om be seen that with the help of Bqs. (2,87) aad 
2.89) y the solutions discussed above can be easily extended 
or batch grinding model *B* also. 

With a view to obtaining sinpler analytical solutions, 
pur j 65-67] has presented a rational approach for a 


I 

I 
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systematic reduction of tlie grinding equation Sq, (2.79), 
wMcla in the limit results in the well known Rosin-Rammler 
type grinding kinetics. E^ur j65j has shown that Bq, (2.79) 
in the cumulative retained mode can be written as 


da.(t) i-i 

“3^ = - B,(t) + RjCt) 

- S. B. .1 

3 1,0 -» 


(2.106) 

where 

i 

Ri(t) = Z M.(t) 
j=1 ^ 

(2.107) 

and 


M^(t) = I^(t) - i^^^(t) 

(2.108) 

For the case when 



(2.109) 


where is some function of size index i only |58 , 77| , 
Eq, (2.106) reduces to 


dil.(t) 

- ' - if— = - a^(t) (2.110) 

whose solution is 

E^(t) - Bji^(O) exp (- Sji^ t) (2.111) 

The condition under which Eq. ( 2 , 109 ) is valid shall be 
discussed later in CSa^ter 5. It may be noted that when 
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cc 

= A , Eq, (2,111) reduces to tlie Rosiu-Ramnl er size 
distribution equation. 

Of several approximate solutions presented by 
Kapur j 65- 67 i » mention should be made of at least one of 
tbem, which is quite in^jortant from the point of view of 
identification of the grinding system. Dividing Eq. (2,106) 
by on both the sides and using two successive interations 

in Cauchy - Picard method of integration of differential 
equations j77 j » Kapur |56) obtained 


R.(t) 
R^ " 


i-1 

s.t + s 
^ 3=1 




where 



r esqp t-1 ^ 

Bj(0) 1 

\ / 



( 2 . 112 ) 




i-1 

S 

3=1 



1,3+1 



(2.115) 


Purther discussion on the utility of this- solution will be 
given in the next section . 

It will be shown in Gh^ter 5 that in general $ the 
size discretized parameters are functions of the particle 
size distribution within the discrete size intervals and 
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therefore, in the strict sense, above mentioned solutions, 
which are based on time - independent parameters, are not 
valid * Only in those cases, where over the time period of 
interest all the parameters can be assumed to be reasonably 
constant, these solutions can find any use • 

2,4.3 Estimation of Time Invariant Size Discretized 

Parameters 

The parameter estimation techniques currently being 
used by various schools can be broadly oLassified as - (i) 
direct determination from the experimental data , (ii) analog 
simulation of the grinding system and cuirve fitting to the 
exp elemental data, (iii) minimization of an error function 
formulated in terms of experimental mass fraction data and 
the predicted results obtained using the analytical solution 
to the size discretized batch grinding eq.uation , and (iv) 
back calculation using the experimental data and the 
analytical solution to themodel eqi nation. 

The most pcpular method for direct determination of 
the rate parameters is to grind single size feeds and plot 
In M^(t) vs. t for the feed size j . Should the resulting 
curve come out to be a straight line, the slope of the line 
gives the value of the rate parameter , If the same 
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experiment is carried out for a small duration of time, 
tlie breakage parameters B. . are obtained by forming the 

j 

ratios of the mass of material that appears in any lower 
size i to the total mass that is ground out of the feed 
size interval 3 , It should be pointed out that accurate 
determinations of bredcage parameters is not possible when 
this method is employed. This is because in an ideal case 
the experiment should be performed for an infinitesimal time 
period such that no particle gets a chance for rebreafcage , 
In practice, the duration of experiment has to be finite 
and, moreoTer, it should be sufficiently large so that 
measurable quantities of material are formed in the lower 
size intervals . G-ardner and Sukan^'najtee j62| proposed that 
10 per-cent of the material may be allowed to bre^ out of 
the feed size interv^ , Hinpel and Austin j6lj and 
Moreira et al* |64| used 1 minute grinding time which 
corresponded to approximately 20 pe2>-cent reduction in 
the feed size materia • Obviously , the parameter estimates 
obtained by these authors are very app2X)Ximate , Austin 
and Ludcie j78( have suggested that several experiments 
may be performed for different grinding times and the 
results so obtained may be ertrspolated to zero time « 

They, however, admit that the procedure is very tedious 
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gfid time consuming , and therefore, it is not ^uoh premising . 

inother variation of this method is to tagg the 
material in the size integral of interest with a tracer, 
which is commonly a radioactive tracer j25, 62, 63, 64( • This 
method has the advantage that the parameter estimates so 
obtained are representative of the actual conditions which 
we desire to simulate . In other words, if the size distri- 
bution of the environment has any significant effect on the 
grinding behaviour of the particles, this effect is also 
taken into account in this procedure * However , since the 
starting feed is not a single size feed but a feed having 
a natural size distribution of interest, the proportion of 
the tagged material may be quite small in some size fractions* 
Hence, the uncertainties in the measurement of broken 
tagged material which reports to the fine size intervals, 
would be relatively large and correspondin^y the estimates 
of the breakage parameters coiig)aratively less accurate • 

Ireeh et al . (52j and Reid and Stewart \5l\ used 
analogue , curve - fitting technique to determine the 
‘best fit* to the experimental data . The set of n 
simultaneous differential equations given by Eq, (2,79) were 
prograjmed for solution on an anaGLogue conputer and the 
values of the parameters were deteimined by trial and error 
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adDustment of the relevant potentiometer knobs. To fit 

% 

the data for n-th size interval (where n is usually 10-12), 
it Would reqLnire to adjust n potentiometer knobs simultane- 
ously , Obviously, for n large (say n > 6) this procedure 
is not practically feasible . Perh^s , for this reason 
, in both the studies some simplifying assui35)tions were 
made to reduce the number of parameters. Preeh et al. 
assumed that the ground material is fed forward to only 
next two or three size intervals. Hence, only three or four 
potentiometers need be adjusted for size intervals 
corresponding to size index i > 4 . Heid and Stewa:^ 
assumed that the breakage parameters are difference - similar 
(noCTalizable) i.e, , B. . = B. . . Hence, only two 

ifj 

potentiometer knobs were required to be adjusted for each 
size interv^ , The sinplifying assuoption made by Preeh 
et al. is obviously quite arbitrary and has no physical 
justification , whatsoever . It will be shown in Oh^ter 7 
that the time invariant difference - similar breakage 
parameter necessarily requires that the rate parameters 
conform to = A • Many exanples can be cited from 
the literature |55,61,71,75,79j which show that the rate 
parameters exhibit a maxima at some intermediate size 
interval , Hence, the assunption of difference - similar 
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breakage parameter is not valid in tbe general case , A lack 
of appreciation of this point has led to inner contradictions 
in many papers | 61,65,64,79 | including the work of Reid and 
Stewart , It can thus be seen that in addition to the 
problem of simultaneous adoustment of a large number of n 
knobs and, therefore, the possibility of introducing a hi^ 
uncertainty in the estimates of those parameters to which the 
solution to the model equation is not quite sensitive, it 
also requires that a large number (at least 2n) / experiments 

A 

should be performed in order that meaningful estimates may be 
obtained in the presence of unavoidable noise in the 
experimental measurements • In conclusion, from the point 
of view of the analysis of grinding systems which do not 
exihibit difference-similar breakage parameters , the method 
of analog simulation is unsuitable and inpractical. 

KlinpeL and Austin j 61 1 and Moreira et al, j64| used 
optimization methods to determine the values of the rate 
parameters • She breakage parameters were obtained using the 
direct method desca:ibed above. The criterion chosen to 
represent the goodness of match was minimization of the 
summed squarred error in values of calculated and experimental 
mass fractions : 

n A 2 

Er •= S [ M. (t) - M. (t) ] 

1=1 . ^ ^ 
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where Er is the error function and M^(t) and M^(t) are 
predicted and experimental oaes ..fractions , respectively. 
Incidentally, Olsen j53a| has also used the same ^proach 
for determination of the selection and distribution parameters. 
Eliiip)el and Austin preferred to use a polynomial in the 
arithmetic mean of ipper and lower sieve sizes of the i-th 
size interval to describe the size variation of the rate 
parameters and deteimined the rate parameters in terms of 
the coefficients of the polynomial . Hiis method has the 
advantage that the parameters can be estimated in terms of 
a fewer number of constants and coefficients. However, in 
order that meaningful parameter estimates may be obtained, 
it is necessary to ascertain that the functional forms chosen 
for the two set of parameters are consistent with their 
basic definition , It will be shown in Ch^ter 6 that the 
time independent rate and breakage parameters are uniqu^y 
interrelated through the functions S(x) and B(v,x) and 
therefore, the functional forms for the two sets of 
parameters cannot be chosen arbitrarily and independently, 

Por exacple, Eliapel aM Austin |6lj assumed the bredcage 
parameters to be difference - similar , In this case , for 
the reasons mentioned earlier , it was then necessary to 
assume - i. instead of a polynomial in , In fact. 
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the rate parameters estimated by Xlinpel and Au^in in 
combination with the assumed difference similar breakage 
parameteis fozm a set, which is not consistent , 

Assuming that the relationship in Eq. (2,109) is valid, 
Austin and luckie |78| showed that the feed size rate 
parameter S-| and the breakage parameters ^ can be obtained 
from the following exact expressions 5 

1 , 

t (2 o114) 


^i,1 


(2,115) 


in [ B, (t)/ R . (0) ] 
in [ B^(t)/ 11,(0) ] 

Herbst and Euerstenau j58j noted that when the 
relationship given in Eq. (2,103) is valid, the rate of 
production of fines is constant for short grinding times 
and in particular 

dE^(t) 


dt 




(2,116) 


where is the cumulative fraction finer than size 

at time t, i, e. 


Pi(t) 


00 

= S M.(t) 
j=-i+1 3 


(2o117) 


m 
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Using the experimentally determined values of and 
and assuming that a , they derived following 

expressions for S. and B. . : 

J d 


S. = 
0 


B. . 


k 


^3 ^ 3+1 

( ) 


a 


(2.118) 


0 


^i 


a 




(2.119) 


■ 3+1 


where kQ and a are system constants determinable from 
measured values . By substituting 3 for i in Eq. (2,109) 
and noting that by definition B. . = 1 , it can be seen that 

dfO 



( 2 , 120 ) 


ilso, it can be readily shown from Bqs. (2,109) and (2,119) 
that 

^i,1 " ^i / ®1 (2.121) 


and 




S- 

3 


( 2 . 122 ) 


It can thus be seen that Eqs. (2.120) and (2,122) provide 
a much more convenient way for obtaining the parameters > 
and B- 4 than what has been suggested by Herbst and 
Puerstenasi • In fact unnecessary oonplications in their 


# 
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work arose because they failed to realize the identity in 
Eq. (2.120) . 

It will be shown in Chapter 7 that the condition in 
Eq. (2.109) does not necessarily inply that a , 

In view of this fact, Bqs. (2.120) and (2,122) are also 
more general than those suggested by Herbst and Euerstenau, 

In addition to the above two criticisms, it should 
be pointed out that the condition S. B. . = Q necessarily 
inplies Rosin - Rammler type grinding kinetics , Erom 
Eq, (2,111) it follows that In Rj_(t) vs. t plots should be 
straight lines in this case, The data for which the ^ove 
method has been applied by Herbst and Euerstenau does not 
exhibit this behaviour as can be seen from ELgure 7 in 
their p^er j58( , Consequently, their estimation of the 
parameters is invalid, 

Justin and luckie j78j have presented a back-calculation 
method, which they have applied for determination of the 
nonr-diff erence-similar (non-no rmaliz able) breakage 
parameters [75 , 80{ . They made two assun5)tions - (i) for 
short grinding periods, amount broken to less than size 
i from breakage of material in size intemral 3 osn be 
approximated by B. . t [ M (0) + M.(t) ] / 2 and (ii) 
while determining B^^^ corresponding to a fixed 3, Bj;^j^(k>3) 
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can be approximated by B. . , . i, e. , the assunption of a 

J— il j J 

difference - similar breakage parameter is approximately 
valid in this case . Suiprisin^y, even though the method 
was proposed for determination of non-difference- similar 
breakage parameters, for testing the method they used 
exact numerically, generated data based on only difference 
-similar breakage parameters (80j . Gonseciuently, it is 
doubtful if the mebhdd will serve the purpose for which it 
has actually ■ been proposed. Moreover, it should be pointed 
out .that the back- calculation methods are extremely sensitive 

m * 

to statistic^ fluctuations in experimental data |49,62l, 
Gardner and Srkanjnaitee j62( have shown using synthetic 
data with controlled random error, that even if the exact 
estimates of breakage parameters are available, the error 
in the estimated values of rate parameters for lower size 
intervals could be as hi^ as 25~50 per-cent , 

Austin and Bhatia |68j and Austin and luckie \ 15 \ 
have proposed some aapirical expressions for parameters 
and B. . which involve five unknown constants. They have also 
described approximate gr^hical procedures to determine these 
constants, which are not very satisfactory. 

finally, mention may also be made of the approximate 
identification procedure suggested by E^ur j 66 - 67 j , fhe 


i 


! 
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method has the advantage that the number of parameters is 
reduced from n (n+1)/2 to only n j66| or 2 n \6l\ , Moreover y 
two grinding experiments with distributed feed are adequate 
for identification of the grinding system . As pointed out 
earlier, the size discretized parameters are inherently 
time-dependent. Hence, if time invariant approximate values 
only are to be used, E^ur*s model should be adequate for 
at least routine simulation purposes , 

We have thus seen that all work in the area of 
parameter estimation has been based on simplified nKDdels, 

In addition to the several restrictive assunptions pointed 
out earlier, parameters have invariably been assumed to be 
time invariant , In Chapters 8 and 9> a novel approach will 
be presented which does not call for making any arbitrary 
assuDptions regarding the nature and structure of the 
parameters . Moreover, it also allows us to calculate time 
variation of the parameters over discrete inteirvals of time, 

2,5 Some Special Grinding Phenomena 

Based on the analysis of experimental data, 
several special grinding phenomena have been reported ^from 
time to time . In the previous sections , mention has already 
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been made of some of these - (1) difference - similar 
breakage parameters, (ii) Rosin - Rammler Kinetics, (iii) 
Alyavdin Kinetics and (iv) a pQwei'.law for size variation 
of the rate parameters . We present below a brief review 
on two of the more important topics - (i) first order 
disappearance kinetics and (ii) zero order production of 
fines f which was not discussed earlier . 

In a recent review, Austin \28\ says that ’ *the 
first order disappearance of material in a discrete size 
interval is a hypothesis which has been proved as a reasoziable 
first approximation It is true that a large number of 
studies can be cited in support of this hypothesis j 25, 48, 
50, 81, 82, 83 I , but at the same time evidences can also 
be produced which go against it [ 41, 42, 84 i . Austin and 
Luckie [ 72 j have shown that this hypothesis is valid for 
a suitably small size interval only . Harris | 45 j has 
also pointed out that the InM^(t) vs. t plots are straight 
lines for some systems and not in case of others , Both 
Harris aud Austin have attributed these deviations from 
the first order kinetics to many factors . Some of these 
are - (i) width of discrete size intervals j 45, 72 [, (ii) 
ratio of the particle size to the mill diameter aud ball 
diameter I- 45, 46, 80 |,,(iii) interactioi^ between the 
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constant rate ,j 31, 42, 58, 83, 85, 86 j, i.e. 

= ^(x) (2.124) 

where ^(x) is the formation rate constant for the material 
finer than size x . This phenomenon was first noted by 
Fahrenwald j 85 j . For a batch ball milling of quartz , 

Arbi tier and Bhrany j 83 ] found that 

<^(x) = (-|t (2.125) 

where x ' and a are size modulus and distribution modulus of 
that portion of the product size distribution which is 
described by the Ga,udin - Schumann distribution, and k^ is 
a constant . Fuerstenau and Somasundaran j 86 | , and Herbst 
and Fuerstenau | 58 j faund that for limestone and dolomite 
alsOythe relation in Eq, (2,125)‘ is valid. 

It has been pointed out that for Eq, (2.124) to hold 
, the grinding of the fine particles t'fcieniselves must be 
negligible | 42 , 58 , 86 | . Mika et al. | 42 } have rightly 
argued that an additional requirement is that the fine size 
fractions do not influence the breakage of coarser material 
in the appropriate time interval. 

Harri s j 31 I showed that for small values. of 

X and t , Eqs. (2.124) and (2.125) can be drived from 

Eq. (2.54) if g(x) is assumed to be of the form 

a 

g (x ) = exp (-x/x*) 


( 2 . 126 ) 
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Mika et al. | 42 ] and Herbst ard Fuerstenau j 58 [ proved 
that the condition given in Eq. (2.109) is sufficient to 
ensure the relationship given in Bq. (2.116), which is the 
size discretized version of Eq. (2.124) . It will be seen 
that the condition in Eq. (2.126) necessarily implies that 
Rosin-Rannnler type kinetics prevails . It ivill be shown in 
Chapter 7 that the condition in Eq. (2,109) also necessarily 
implies Rosin-Rammler type kinetics , It then follows that 
the systems exhibiting the phenomenon of zero-order production 
of fines must conform to Rosin-Rammler type of kinetics, 

which surprisingly is not the case in both '■ studies 

. , I 

-§£-4 Mika et al.^and Herbst and Fuerstenau 4 *r 6 ’ . 

In a later publication j 45 I , Harris pointed out 
that the conventional method of plotting F(x,t) vs. t on an 
arithmetic grind is insensitive, and a close inspection of 
the grinding data reveals that it exhibits a systmatic 
departure from Eq. (2.124) .He suggested that the following 
equation seems to be more appropriate for describing the 
data 

F(x,t) cc x“ (2.127) 

A formal derivation of Bq. (2.127) from Bq. (2.65) has been 
givenby Harri s in an earlier publication | 38 | . Employing 
data from various sources, Harris has shown that the constant 
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*a* normally lies in the range - 0.5 < a < 0.3 . 

It is quite evident from the above discussion that 
the two hypotheses of first order disappearance of a discrete 
size fraction and zero-order production of fines, are valid 
only in a very approximate sense . Harris | 45 1 is right 
in proposing that these should be discarded as the general 
hypotheses . Further discussion on these two phenomena shall 
be given in Chapter 7. It may be mentioned briefly that the 
observed deviations are only natural consequences of violation 
of the Rosin-Rammler tsrpe of grinding kinetics . 

2.6 Empirical Correlations for Variation in Rate Parameters 
With Various Operating Variables 

A knowledge of the quantitative correlations for 
variation in the model parameters is quite important from the 
point of view of simiilation studies and also for establishing 
the optimal working conditions. lo date the effects of 
particle load j 25, 42, 82, 87 - 91 I’ , ball load j 91-93 |, 
mill speed j 92 [ , ball diameter j 82, 94 - 98 |. , grinding 
media shape { 99 j and grinding media density j 93 | have 
been investigated . A brief account of these findings is 
presented below . 
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Bowdish I 82 j found that the rate parameter S decreases 
as the hold-up of the particulate mass W is increased r. Bor W 
corresponding to 100 per-cent filling of the interestitial 
volume of the static ball charge, the value of the absolute 
grinding rate SW v^as f)und to be maximum . This result is in 
agreement with the findings of Kapur aod Agrawal { 88 ] 
and Mika et al. j 42 [ . ILjman | 89, 90 U ^^^i'ka et al. \ 4-2 j 
Kelsall et al. j 56, 87 j and Austin et al. j 91 | have shown 
that the rate parameter was inversely proportional to the 
particulate mass W when the particle loadings exceeded 100 
per-cent filling of the interstial volume of the static ball 
charge , In the unsaturated regime (that is, less than 100 
per-cent filling ), the rate parameter was found to be 
approximately inversely proportional to the square root of 
•the particulate mass | 42 | . However, a single correlation, 
valid for all size intervals andfor both starved and saturated 
mills, has not been presented so far . In particular, it 
has not been explicitly demonstrated whether the particulate 
mass affects the rate parameters for all size intervals of 
interest to the same extent , That this is indeed the case 
we shall show in Chapter 10 . 

Herbst and Puerstenau ] 92 j found the rate parameter 
to be independent of the ball load, when the latter was less 
than 40 per-cent of the load required to fill the mill 
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completely and the particle load corresponded to 100 
per-cent filling of the interstices of the ball mass at 
rest , For higher ball loads the value of the rate 
parameter was found to decrease, which was attributed to 
ball-ball interactions . The results of Herbst and 
Fuerstenau ] 92 j as well as Austin et al. j 91 1 show 
that the absolute grinding rate SW increased approximately 
linearly with ball load for small loads and attained a 
maximum at a load corresponding to 50 per-cent of that 
required for complete filling of the mill . Increasing the 
ball load beyond this percentage did not result in any 
change in the absolute grinding rate . The broad trend of 
results obtained by Kelsall et al. 1 93 | is also the 
same , 

Herbst and Fuerstenau ] 92 | have reported 
results on variation of rate parameter with mill speed 
at constant particle and ball loads , Rate parameter was 
found to increase approximately linearly with increasing 
speed of rotation, attaining a maximum at about 80 per- 
cent of the critical speed, and then decreasing sharply 
as the mill speed approached 100 per-cent of the critical 
speed. According to Rose { 100 | , at slow speeds of 
rotation the dynamic configuration of the ball mass 
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changes little with increasing speed and therefore, the 
rate parameter on a per revoluation basis is expected to 
remain essentially constant for these conditions as has 
been observed by Herbst ani Pnerstenau . Ba.sed on a mathem 
-atical theory of ball motion due to Davis j 101 | , for 
the experimental conditions employed by these authors* 
the best operating speed is calculated as 80 percent of 
the critical speed. This also agrees well with the 
experimental findings of Herbst ajid Puerstenau . 

Bowdish j 82 j studied the effect of ball size, 
keeping total weight of balls constant, on the rate 
parameters for different sizes of limestone feed . He 
found that the rate parameter were directly proportional 
to the surface area of the balls . However, later studies 
I 89, 90, 94,-98 ( have confirmed that for a given mill, 
material, particle and ball loads, there is an optimal 
ball size for the maximum i*ate of grinding of a given 
type and size of particles. That is, for a given particle 
size, the rate parameter first increases with increase 
in ball diameter, reaches a sharp maximum and then 
decreases « Moreover, the curves relating rate parameter 
to ball diameter have been found to be similar in shape 
for all particle sizes, and both the optimal ball size 
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and the corresponding value of the rate parameter have 
been found tj increase with increasing particle size. 

The portion of the curve in which the rate parameter 
decreases with increase in ball diameter has been found 
to be of the foira j 96 [ and -’C i 94 1 , 

where B is the ball diameter,. However, no general 
expressions are available relating the rate parameter 
as function of ball aixL particle size. 

In addition to the studies summarized above, 
Kelsall et al, have reported results on the effect of 
grinding media shape j 99 j and grinding media density 
1 93 1 also, Hovyever, the data presented is not sufficient 
for deriving quantitative correlations. 

It should be mentioned that of all the variables 
discussed above only mill speed ani grinding media shape 
have been reported to have any significant effect on 
the breakage parameters j 92, 99 | . Again, no quantitative 
correlations are available for describing the variation 
of bresikage parameters with these two variables. 



CHAPTER 3 


OUTLINE OE WORK 


The present work can be conreniently described 
under seven subtitles as follows : 

(I) 4 Pseudo- Similarity Solution to the Integro- 
Differential Equation of Batch G-rinding 

With a view to obtaining a convenient closed 
form analytical solution to the fundamental integro- 
differential equation of batch grinding for more general 
forms of the selection and breakage functions than those 
used hitherto by a number of invest i gat iors j 14, 24 , 30 , 
31 I , a new solution has been obtained using Kapur's sim- 
ilarity solution and standard mathematical transformation 
techniques . The resiilting size spectra esiiibit a non- 
self-preserving nature, which has been illustrated by 
simulation studies . The solution is the most general 
close form solution available at present „ 

(II) Exact Mathematical Definitions of the Parameters 
in the Discrete Size Grinding Models 

As pointed out earlier in Chapter 2, most of 
the investigators have either formulated the discrete 
size models without any reference to the more fundamental 
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size continuous model and the two basic functions S(x) and 
B(v,x), or they hare failed to obtain the correct relationships 
for the discrete model parameters in terms of the two basic 
functions « This has led to considerable confusion and 
misinterpretation of many experimental data . With a view 
to providing a proper framevfork for formulation of the discrete 
size models, the exact mat heme" ical definitions of the 
discrete size parameters have been derived from the fundamental 
integro-differential equation of batch grinding using a formal 
mathematical procedure . It has been shown that all the 
parameters are interrela.ted in a unique manner through the 
functions S(x), B(v,x) and M(x,t), and hence, cannot be 
assigned any arbitrary values, independent of each other . 

Anew function L(v,x) = S(v) B(v,x), which 
plays a central role in the analysis of discrete size models, 
has been identified and its important properties have been 
delineated. 

(Ill) Realizations of Various Specialized Discrete Size 
Grinding Models 

It is possible to simplify the expressions for 
size discretized parameters by specializing the functional 
forms of the functions M(x,t) and I(v,x). Three such 
specialized cases have been identified and the important 



89 


oliarac'teris’tics of the parameters of the resulting models 
have been delineated » 

(lY) Analysis of Some Special G-rinding Phenomena 

The phenomena of first order disappearance kinetics 
, zero-order production of fines, Alyavdin kinetics, difference 
similar breakage parameters, and variation of the rate 
parameters with particle size according to a power law, 
have been investigated in light of the structure of the three 
specialized discrete size models proposed. The pertinent 
discrete size models associated with these phenomena have been 
identified and delimited in this context, 

(Y) Identification of the Grinding System in Terms of a Mill 
Punct ion 

A novel mill function approach is presented for 
identification of a linear ani time- invariant batch grinding 
system, in which size reduction occurs by impact fracture, 
chipping and abrasion. The approach is based on a size 
continuous and discrete time representation of the grinding 
system. The proposed identification scheme ha,s been demonstra- 
ted for a batch ball mill - dolomite grinding data , 

(YI) Estimation of Size Discretized Parameters 

A ‘1 ,0* variable suppression technique has been 
devised for exact back calculation of all the size discretized 
parameters S* and B. j from a kn.owledge of the input-output 
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matrix 0 in (2.103), Using the estimated mill functions 
for two sets of batch grinding data for dolomite, the 0 matrices 
have been generated over several successive intervals of 
fixed grinding time and, for first time, the time variation 
of the model parameters due to their dq)endence on the size 
distribution of the particles within the individual discrete 
size intervals, has been successfully confuted. 

Based on mill-function a new set of parameters have 
be@3 defined and named as ’reference parameters' • inhese are 
useful from the point of view of conparing different grinding 
systons. 

(vii) Empirical correlations for Tailation of Sizo - 

Discretized Bate Parameters With Particulate Mass 
and Ball Size 

Systematic gr^hical procedures for determination 
of the enpirical corrriations for the variations in time 
ind^endent rate paraineters (Bass Model ) with particulate 
nass and the size of the grinding m6dia ball have been 
ieveloped and successfully demonstrated for two actual 
grinding systems. 
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OHiPlER 4 

A PSSUroO-SIMIIARITI SOIUiDION TO THE IITEGRO-HIFEEREiVTIAIf 
EQUATIOI OE HATCH GRIHHIHG 

We have seen in Chapter 2.3.3 that it has not been 
possible to obtain a convenient closed form analytical solution 
to the fundamental integro-differential equation of batch 
grinding (Eq. 2.36) for general forms of the selection aM 
breakage functions . On the other hand , the specialized 
mathamatical forms of the selection and breakage functions 
used by Gaudin and Meloy j24| , Austin et al. |50j , Harris 
|31| and Kapur j 14] seem to be much too restricted to be 
realistic for a given mill-nBt erial systan. In this Chapter 
.an attempt has been made to extract solutions to Bq. (2.36) 
for a more general class of the selection and breakage 

functions . The strategy involves a suitable mathanatical 

eu _ 

transformation of particle size x into a psuedo-ssize u, a 
similarity solution to the grinding equation in u , followed 
by back transformation into x, resulting in a pscu^o-similarity 
solution to the grinding equation in particle size . 

4.1 Psaudo- Similarity Solution 

Let the grinding functions be given by 

SCx) = A [ g(x) ] 


(4,1) 



and 


(+- 2 ) 

wiLere g ig any moilotonically non- decreasing function in 
particle size such, that g(0) = 0 . Substitution of Eq^s. (4.1) 
and (4.2) into Eq. (2.36) gives 


aM(.| ^ t2 = - ^ [ g(x) M(x,t) 

+ S A[g(v) f [ f[|} M(v,t) dv 
° (4.3) 

We seek a solution to this equation by a suitable one to 
one transformation of the variable x „ let us define a 
transformed variable , the pseudo- size u , in the following 
manner 

u = g (x) (4.4) 

and w = g (v) (4. 5) 

Hence , the distribution in psx'.ado-size I(u , t) du , 
defined as the mass fraction of particles with attribute in 
the range of u to u + du , is related to M(x,t) dx by 


M(x,t) = Y(u , t) I ^ j,. 


or 

M(x,t) = T [ g(x), t] ( I 

Moreover 



1 )^ 
8u w 


dg(x) 

dx 


(4.6) 

(4.7) 


(4.8) 
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Substitution of Eqs. (4.4) - ( 4 . 8 ) into Bq. (4.3) gives 
= - A Y(u,t) + / A / L ( I Y(5 ,t) d 


u 


w 

(4.9) 


which has a similarity solution 1 14 




(4. 10) 

where 

CXJ 

2 ^^(t) = / u Y( u,t) du ( 4 . 11 ) 

' 0 

Substitution of Eqs. (4.4) and ( 4 , 6 ) gives 


= / g(x) M(x,t) dx (4.12) 

Using Eqs. (4.4) and (4.7) , we back transform the distribution 
in u into the distribution in particle size x and obtain 


M(x,t ) 



C- ffi" ( 


)“i 


^g(x) 

dx 


(4.13) 


which is the desired solution to Eq, (4.3) . It is shown in 
the Appendix 2 that when g(x) 5 = x^ , we obtain the original 
similarity solution of Eapur from Eq. (4.13) , thereby 
confirming the validity of our approach , It should be noted 
that while the distribution in u , Y(u,t) , is self-preserving, 
the distribution M(x,t) in scaled variables x/x/^(t) or g(x)/ 2 ^(t) 
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is no longer self-similar because of the dg(x)/dx term present 
in Eq. (4.13) . This term is independent of x onlj when g(x) 
is a linear function of x , Thus under a linear transformation 
in Eq. (4.4) , the size distribution in g(x)/2^'^(t) retains 
its self- similar character . 

4.2 Simulation Results 


In order to illustrate the shape of the computed 
size distribution curves of M(x,t) , we have a.ssumed the 
following two forms for the function g(x) 



g^(x) = 0.5 [ X + 10~^ x^ ] 

(4. 14) 

Hence 

dg.(x) , 

= 0.5+ 10 X 

(4.15) 

And 

gg(x) =10^ [ ezp (10"^ x) - 1 3 

(4.16) 

Hence 

dgpCx) ^ _3 

^ = exp (10 ^ X ) 

(4.17) 

where particle size x is expressed in miciwns , 

ClteKUlative 


size distributions were computed in u domain using Iq. (4.10) 
with A = 1, a = 0.5 , p = 1, and 7 /^( 0 ) = 5657 , for four 
times corresponding to 2/^ =2>'^(0) , ^(0)/4 , 2->j(0)/l6 , 

and 2^-|(0)/64 .. Distributions in x domain were obtained by 
transforming u into x using Eqs. (4.4) » (4.14) and (4.16) , 
Figures 4,1 and 4.2 show these distributions at various times 
on log-log grid . The shapes of the distributions are quite 
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similar to those based on numeirous experimental date reported 
in the literature * 

In order to illustrate the nature of deviation from 
the self -preserving behaviour , we recall jl4j that the shape 
of the similarity distribution in u domain is time invariant 
when the cumulative distribution function is plotted as a 
function of u/uq ^ , where Uq ^ is the median size , Hgures 
( 4 . 3 ) and (4.4) for transformations g^ and » respectively, 
shovi/ that the curves in particle size x do not collapse and 
are significantly non-similar . The dotted curves are the 
self - preserving forms in psott^o-size u domain , The nature 
and the extent of deviation from a self-similar shape depends 
on the type of transformation ertployed . With increasing 
grinding time more and more material concentre.t es in the finer 
size range , Inspection of Eqs. (4.14) and (4.16) shows that 
as X approaches zero , the transformations tend to become 
linear , Therefore , the distributions in x should tend to 
the self-preserving form T\/ith increasing grinding time . This 
is readily confirmed from Eigs. (4.3) and (4.4) , where it 
can be seen that pstudo- similarity curves eventually merge 
with the self-preserving curves in u domain . 

,3 Discussion 


The frequency function of the cominuted particles in 
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Eq. ( 4 . 13 ) is the most general analytical expression ever 
derived from a formal kinetic grinding model. It is much more 
elastic and versatile than the Ro si nr- Raninl er and the generalized 
gai-fflua distributions . The equation describes the complete 
trajectory of the particle spectra in time , apart from a 
small initial period , As such its inportance , specially for 
simulation studies , is self evident . The combined transfor- 
mation - similarity solution method presented here is of 
great utility | it can be employed for obtaining new solutions 

to a variety of particle population problems 1 103, I in the 

104 

way similar to the one described above . . 
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GHiiPTER 5 


EX4CT MAIhBdjiailGilj DEFINITIONS Oi' THE PiilUviETERS 
IN THE DI3GHETE SIZE GRINDmG MODELS 


In view of tHe discussion in Chapter 2 on the current 
status of the literature on discrete size grinding models , 
it is quite evident that the establishment of the correct 
mathematical definitions of the discrete size parameters in 
tems of the two basic size continuous f unctions j,S(x) and 
B(v, x). j is a task of some importance , In this chapter 
therefore , the required expressions for the parameters of 
the two size discretized models have been obtained using a 
fomal mathamat ical procedure , which consists of first an 
appropriate size discretization of the fundamental integro- 
differential equation of batch grinding and then a team 
by team comparison with the two master kinetic equations for 
discrete size models . 

5.1 Derivation of the Exact Mathematical Definitions 
of the Parameters for Model 'A' 

Substitution of Bq. (2.34) iu the following 


identity 
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dt 


d [ - I.(t) ] 

dt 


(5.1) 


gives 

dMj;(t) 

“”dt 


"0 


/ S(v) B (v, M(v, t) dv 


^i-1 


X, 


0 


- f S(v) B (v 5 M (v, t) dv 


(5.2) 




Break up of tlie integrals- over discrete size intervals and 
rearrangement of the resulting terns gives 


dM^(t) 


^i-1 


dt 


= -/ S(v) B (v, x^) M(v,t) 


dv 


^i 


i-1 ^0-1 


+ 2 / S(v) f B(v,x. J-B(v,Xj_)]M(v, t) dv 

0=1 X. 

3 (5.3) 


Term by term comparison of Eq.s. (2. 79) and (5.3) shows that 


X 


i-1 


S^(t) M^(t) = / S(v) B(v,xp M(v,t) dv 


^i 


Hence , the instantaneous rate parameter is 


(5.4) 


5i(t) = 5(v,Xj.) M(v,t) dv (5.5) 


^i 


We note that the mass fraction of solids in size interval i is 
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^i-1 

= / M(x,t) dz , 1=1,2, ... (5.6) 


Hence , an alternative equivalent e3q)ression for the rate 
parameter in teniis of the size continuous functions can he 
wirtten by combining Eqs. (5.5) and (5.6) 

^i-1 

/ S(v) B(v, X .) M(v,t) dv 

S (t) = r (5.7) 

^ ^i-1 

f M(v, t) dv 


In same manner , the expression for the instantaneous 
breakage parameter is 

/ S(v) [ B(v, X. J - B(v,x.) ] M(v,t) dv 

X. 

_J 

/ S(v) B(v, X.) M(v, t) dv 

X ^ 

3 

, 3 7^ i (5.8) 




Substitution of Eq. (5.8) in Bq. (2.82) gives the following 
expression for the instantaneous cumulative breakage 
parameter 






s S(v) B(v,Xj_) M(v,t) dv 




/ S(v) B(v,x.) M(v,t) dv 

X. ^ 


(5.9) 
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The relationships in Bis. (5.7) , (5.8) and (5.9) are of great 

importance , Much of the controversy and confusion - frequently 

leading to erroneous conclusions - prevailing in the recent 

grinding literature may be attributed to our failure to 

recognize these relationships and derive the correct 

conclusions therefrom , It should be noted that Bis. (5,7) 

and (5.8) clearly establish that - (i) the parameters 

and B. . are not independent of each other as both are 
-L? J ^ 

determined by the three functions S(x) , B(v, x) andM(x, t) 
jointly , (ii) there is no one - to - one physical 
correspondence between S(x) and or between B(v, x) and 
B. . , and (iii) the parameters are functions of the 
instantaneous particle size distribution and , hence , vary 
with time , In view of the last abservation , it becomes 
necessary to rewrite the Bi. (2.79) as follows 


dM^(t) 

dt 


i-1 

S.(t) M.(t) + 2 

^ ^ j=1 


S^(t) B.^ .(t) M^.(t) 


, i = 2 , 3, (5.10) 


It is interesting to note that even though the grinding 
systen has been assumed to be linear and time - invariant , 
the size discretized representation of the kinetic model 
is that of a time varying system , This apparent anamoly 
can , however , be easily resolved , It should be recalled 
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that the model parameters refer to a size range of particles 
and not a sin^e particle size , Also, the contribution to 
the value of a model parameter of each particle size within 
a discrete size interval , is proportional to the relative 
mass density of particle population at that size . Since, 
the latter varies with grinding time , it is easily seen 
that the time variability of the model parameters is 
inherent in the discretization procedure itself and , 
therefore , cannot be avoided , 


I 


5.2 Derivation of the Exact Mathanatical Definitions 
of the Parameters for Model 'B' 

Let us discretize the integral tem on the 
right hand side in Eq. (2.36) as follows 

I 

c'x 

-X. 

1 ^D-1 


= - S(x) M(x,t) + / 's(v) -®^-^-M(v,t) dv 


+ 2 s' S(v) M(v,t) 

j=i-1 X. 


) dv 


(5,11) 


Now, we integrate Bq. (5,11) on both sides over a size 
interval 2 :j_ < x < x^__^ , hence 


^i-1 

8 f M(x, t) dx 
^i 

8 1 


^i-1 

/ S(x) M’'(x,t) dx + 
^i' 
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^i-1 




•S(v) M(v, t) dv dx 


^i-1 X. 


■^3=Li 4 I ’ ^ 




(5.12) 


Now, using Eq, (5.6) and carrying out a tern by tem 
comparison with tiie model equation Eq. (2. 86 ) » we obtain 


1-1 


s^ (t ) = 


/ S(x) M(x,t) dx 


/ M(x,t) dx 


(5.13) 


^i-1 ^i-1 


/ S(v) 


gi n ? . 

9x 


M(T,t) dv dx 


.(t) 


^i-1 


/ S(v) M(v,t) dv 


(5.14) 


^i-1 

/ / S(v) 


M(v,t) dv dx 




^ 3-1 

/ S(v) M(v,t) dv 


,i 9^ 3 (5.15a) 
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wliioh cazi be simplified as 




/ sir) [ B(v,x^_^) - B(v,x^) J M(v,t) dv 
^D 


/ S(v) MCvjt) dT 


, i 7^ D (5.15) 


We recall that by definition 


b. . (t) = 1 - E b. (t) 
0=i+1 


(5.16) 


Thus , by substituting Eq. (5-15) in Eq.. (5.16) , a relatively 
single expression is obtained for the inner breakage parameter 

^i-1 



(t)=1- 


/ sir) B(v,x^) M(v,t) dv 


(5.17) 


/ S(v) M(v, t) dv 


It can be confinned by carrying out some mathematical 
manipulations that Bqs. (5.14) and (5.17) are identical . 

A comparison of the esqpressions derived above with 
those mentioned by Whiten j54( will reveal that his 
equations are grossly incorrect . Similarly, by substituting 
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the relationshiiB used by Olsen j53j for S(v) and -BCvjx) 
and evaluating the integrals in the expressions for s^^, 

^i,i ^i,j ’ shown that these results do not 

match with the corre^onding expressions derived by 01senj53j 
and Used by Olsen and Krogh [Tlj . Cfonsequently , the 
analysis of the grinding operation as presented by these 
authors is suspect ani needs a revision . 

It may also be anphasized that both the rate and 
breakage parameters S. and 3. . , and , the selection and 

1 J 

distribution parameters s. and b. . are time dependent i.e, , 

1 J 

the latter set of parameters are also functions of the 
particle size distribution within the individual size 
intervals . 

Henceforth in this work we shall restrict ourselves 
mainly to the more popular model i.e., the model 'A’ , 

5.3 Ihe Selekage Itmction 

¥e note that in the expressions derived above for 
model *A’ , the selection and breakage functions occur in 
association only , The product term S(v) B(v,x) , which will 
appear repeatedly in the discrete size analysis , deserves 
a name of its own • We shall call it as the ''selekage 
function X(v, x)^ that is 

L(v,x) = S(v) B(v,x) (5.18) 
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We mention in passing that the differential fom of the 
selekage function 

= S(v) 5(v,x) (5.19) 

= G(v,x) 


first appeared in the integro - differential equation of 
hatch grinding formulated hy Bass )22j . 

We now present some of the more important properties 
of the selekage function , Conservation of mass requires that, 
the breakage distribution function be a nomalized one , 
hence 

B(v,v) = 1 (5.20) 

Itom definition of the seldcage function in Eq. (5.18) we 
have . 

S(v) = L(v, v) (5.21) 


and 

B(v,3:) 



Moreover , 


= ilV 


9Ii(v.x) 

ax 


( 5 . 22 ) 

(5.23) 


In principle then the knowledge of the selekage function 
is equivalent to the knowledge of the selection and 
breakage functions . 
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Tile cumulative breakage function B(vjx) is a 
monotoni cally non-decreasing function of x and if S(v) can 
assume only positive definite values , it follows that the 
selekage function is also restricted to positive definite 
values ,only . Moreover 

B(v, 0) = 0 (5.24) 

therefore 

l<(v,0) = 0 (5.25) 


5,4 The Master Equations in Terms of the Seldcage function 

Tbr the reasons mentioned above , in the analysis 
of the discrete size grinding models , it is preferable to 
write the mathamatical expressions in terms of the selekage 
function , Eor ready reference , it is therefore convenient 
to rewrite some of the important equations replacing the 
product term SB by 1 , Consequently , the forms of Sqs. (5.3), 
(5.7) j (5.8) , (5.9) and (5.17) are , respectively . 


dM. (t) ^i-1 

— Ig — = / I.(v, x^) M(v,t) dv >. 

Xi 

+ S I ^ [L(v,x:^ -j) - L(v,x^)3 M(v,t) dv 

3^=1 X. 

(5. 26) 


I 

[ 


[ 


i 





(5,27) 


/ Ii(v, X. ) M(v,t) dv 

X. ^ 

X 


f M(v, t) dv 


^D-1 

/ [LCv,Xj^_^) - L(v,x^) ] M(v,t) dv 


^ 3-1 

/ L(v, X. ) M(v, t) dv 

■ ,i ^ 2 (5.28) 


Ii(v,x^) M(v, t) dv 

^ (5.29) 

3-1 

S I/(v, X. ) M(v,t) dv 


^i-1 

/ Ii(v,Xj_) M(v,t) dv 


S S(v) M(v,t) dv 


1 


(5.30) 
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GFLaPIER 6 

RMIZA[DIOi\'S OP VARIOUS SPECIALIZED DISCRETE SIZE MODELS 

Eor an inherently linear and time- invariant systsa 
tile general most discrete size model is given by Eqs. (5. 10), 
(5.27) and (5.28) in which, both sets of grinding paramet ers 
are non- unique because of their dependence on the instantaneous 
particle size distribution , It is recognised that the 
grinding equations for the G-eneral model present iormidable 
mathaaiatical complexities . therefore , there exists a 
geniune need to e^q^lore the possibilities of simplification 
in the grinding model . In this chapter three such cases have 
been identified and the characteristics of the resulting 
specialized discrete size models have been delineated , 

6.1 Ihe Bass Model 

The Bass model, which is given by Eqs, (2,40) , 

(2.72) and (2.73) , is characterized with time-invariant 
parameters v. . , In alternate derivation of this model is 
possible using somewhat different arguments . Let us assume 
that we are justified in approximating the particle size 
distribution within the discrete size intervals by a unifom 
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distributioja . This implies 


= d. (t) 


, X. < X < 


( 6 . 1 ) 


dj_(t) 


M^(t) 

^i~1 " ^i 


( 6 . 2 ) 


Substituting these ecLUetions into Eqs. (5.27) and (5,28) , 
the following results are obtained 


/ L(v,x. ) dv 




(6.3) 


S [ L(v,x^_^) - L(v,x^) ] dv 




/ I(v,x . ) dv 

"3 


(6.4) 


Hence , 


®3 ®i,3 = 




^i-1 - ^i 


(6.5) 


Prom Eq. (5.19) it follows that 


Ii(v5X. ) = / S(v) B(v,x) dx 

^ 0 


( 6 . 6 ) 
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Substitution of Eq. (6,6) into Bq. (6,5) gives 


X. 

1 

/ 

0 


\-1 


/ S(v) B(v, x) dv dx 


X. 

1 


^i-1 “ ^i 


(6.7a) 


Discretizing the first integral Eq., (6,7a) can he rev«T:itten 


as 


'^k-l 


k > i Xi, 


/ / S(v) B(t, x) dv dx 




X. 


i-1 


^i 


(6.7) 


following Bass l22| , let 


^-1 ^ i-1 


/ / S(v) B(v,x) dv dx 


k,i 


“ ^i 


( 6 . 8 ) 


We have 


^i ^ ^k i 

1 k > i 


(6.9) 


Similarly we can show that in Eq. (6,5) 


^i-1 


S. B. . 
D 


/ / S(v) B(v,x)(irdx 

^D-1 - ^3 


( 6 . 10 ) 
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Hence 




( 6 . 11 ) 


Substitution of Eqs. (6,9) , (2.72) and (6,11) into the 
master equstion Eq. (5.10) immediately leads to the size 
discretized time invariant parameter grinding model in Eq. (2.73) 
i.e. , the Bass model . 

Bass j22| had suggested that the inequality in 
Eq. (2.74) concurrently implies a vanishin^y small inner 
breakage parameter b. . . In order to check this assertion 

Ij X 

we substitute Eq. (6.1) into Eq. (5.30) 


/ I(v, 2 :. ) dv 

•V ^ 


/ S(v) dv 


( 6 . 12 ) 


and to obtain some numerical values of b. . , let us assign - 

Xy X 

simply for the sake of convenience in ccanputation - the 
widely used forms for the selection and breakage functions 
as given in Eqs. (2.56) and (2,57) . therefore , the 
selekage function is 

l(v,x^) = A v“”P 


(6.13) 
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Substitution of Eqs. (2.56) and (6.13) into Eq. (6.12) gives 



1 - [ 


g + 1 

a - p + 1 


] [ 


6 “+'' -1 ^ 


1 \ ] IBS , a-p. 1 

6 - 1 


, a - p + 1 0 

(6. 14a) 

= 0 (6.14b) 


wiiicia is independent of tbe size index . fable 1 gives 
values of the inner breakage pr.a?«-ivtlr for various eombincit ions 
of g and p for two sieve size ratios 6 = fZ and ^f2 when we 
impose the Bass model on the grinding sjstan . 

It is clear that in general the inner bre^-cage 
parameter is not insignificant and a substantial aciount of 
material on breakage falls back into the original size 
interval even when the sieve size ratio is only ^2 . The 
value of the inner breakage parameter is quite insensitive 
to the variation of the selection function with particle 
size but increases significantly with the slope of the 
breakage function , which is intutively reasonable , It 
can be concluded then that the sieve size ratio 6 as snail 

as ^2 is not adequate to satisfy the Bass condition , 

However , as shown above , in order to generate the 
Bass model it is not necessary to stipulate that b- . = 0 , 

The sufficient condition is that the particles are uniformly 



distributed within a size interval . it io j.-, . 

•i-vaj. . i-c IS true that particles 

cannot ranain mlfomly distributed within the size Internals 

throu^out the course of grinding , however , as a matter 

of practical necessity , we may arbitrarily Impose a piecewise 

uniform distribution . This approximation becomes increasingly 

valid as the size inteival width decreases . ilso it is 

physically more realistic because our approximation does not 

ignore the fact that the daughter particles indeed fall back 

into the original size interval . 

As pointed out earlier in Chapter 2 the Bass 
model in Eq. (2.73) can also be represented by Eq, (2.79) , 
which is the most widely used discrete size grinding equation 
in the literature | 42 , 49, 50, 55 - 67 J . However , two 
important differences should be noted in respect of the 
model fomulation , Eirst , unlike in the works of Mika 
et al. j42| and Herbst and Euerstenau {58j , the derivation 
of Eq. (2.79) does not require the assumption of b. -=0 . 

Second , paraneters and ^ , as given in Eqs. ( 6 , 3 ) 
and (6.4) , are not uncorrelated or independent of each other 
as regarded by most of the workers [ 42 , 49 _ 52 , 55 , 58 , 61, 

62, 73 j . Both sets of paraneters are functions of the 
selekage function and thus, are interrelated in a unique 
manner , It also follows that these parameters cannot be 
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ec[uated with mutually independent phenomena of selection for 
breakage and breakage distribution function as originally 
proposed by i^stein j16j . 


6.2 The RSI* Model 


A inspection of Eqs. (5.27) and (5.28) shows that 
it is not essential to assume a unifoim distribution within 


a size interval , as done in the previous case , in order 
to generate a size discretized model with time invariant 
parameters ,• In some cases if it should turn out that the 
selvage function L(v 5 x) is independent of variable v, that 
is 

l(v,x) = A Lg (x) (6.15) 

where L 2 is seme function of x only , then it can be easily 
shown that the form of this reduced selekage function is 
necessarily as given below 


Ii(v, x) = A 


I'2(^) 


(6.16) 


where now A > 0 is a constant . 

Thus from Eqs. (5.21), (5.22) and (6,16) 


S(v) = 

A l2(v) 

(6.17) 

B(v,x) 

12^^) 

1-2 

(6.18) 
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Since b ( v , x ) is a moncbonically non- decreasing function of 
X , it follows that IgCx) and , hence , S(x) is also a 
monotonically non- decreasing function of x , Again, combining 
Eq. (5.24) with Eq. (6.18) gives 

1-2^0^ " ° (6.19) 

therefore 

s(0) = 0 (6.20) 


Some of the results mentioned above were derived earlier by 
Austin et al. j30j . 

Substitution of the reduced selekage function 
(ESE) into Eqs. (5.27) and (5.28) gives 


Si = A Ig (Xi) 

and 

^ Lg(x._.,) - Lg(x.) 


( 6 , 21 ) 


( 6 . 22 ) 


Substitution of these time- invariant parameters into 
Eq. (5.10), the master equation , leads to the size 
discretized ESE model as follows 


dM^(t) 

dt 


= - Al2(x.) M^(t) + 


C I-2Ui-i) - 1-2 Ui) ] 


i-1 
2 M. 

0=1 3 


(6.25) 
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Ihe inner breakage parameter for the RSP model is 


1-1 


M(v, t) dv 


b._.(t) = 1 - 


1-1 


(6.24) 


f S(v) M(v,t) dv 


which is not time - invariant . It will be seen that if 
a = p , Eqs. (2.56) and (2,57) reduce to RSE model with 
~ » ^^<3. S(x) = A x'^ , Now , Kapur's similarity 

solution |14| can be used to compute the time variation of 
b- . . Substitution of Eqs. (2.59) into Eq, (6.24) gives 

X, X 

an esplicit expression for b. . (t) as follows 


= 1 - 


^1 1 

X- J V 

X. 

1 


ST” 1 


X 


X~1 


a 




(6.25) 


Eor simulation was assigned a value 1 and i, a, 6, and 
were chosen 1, 1 , f2 and f 2/4 , respectively . 

Units of time used were .^minutes . Eor illustration of the 
results , has been denoted by 10 mesh sieve size . It is 
evident that in theory the similarity solution covers the 
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range 0 < x < ~ , therefore , the choice of the top size 
interval in this case has to he little arbitrary . This , 
however , does not affect the accuracy of the ccmputations , 
since the size distribution over the entire size range 
is uniquely described by Eq. (2.59) , and is independent of 
the choice of the top size interval , 

The computed size distribution curves at the 
start of the experiment and at t = 9 minutes are shown in 
Eigure 6,1 , The variation of the inner breakage paraneter 
with grinding time for four size intervals i = 1, 3, 6 and 
10 has been shown in Eigure 6.2 , It can be seen that , at 
least , in case of cogfecse size intervals , the variation is 
quite significant . Moreover , the value of b. . for 
ooQrse size intervals is markedly different then that for 
the same in the Bass model . 

6. 3 liie SSE Model 

In this model we f omul ate a. separable selefcage 
function (SSB) as follows 

L(v,x) = li^(v) (6.26) 

where L^(v) is some function of particle size v, 

It can be shown that in case of the separable 
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selekage function the fom is necessarily 
L(v, x) = h(v) 


l2(x) 


where h(x) is any positive definite function over the 
range 0 < x < x^ . It follows that 


(6„27) 


(v) = 

bCv) 

L^) 

CM 

o 

and from Eqs, (5.21) 

j (5.22) and (6.27) 


S(v) = 

h(v) 


and 



B(v,x) = 


(6.29) 


It should he noted that in contrast to the RSP model j 
in this case S(0) need not necessarily be equal to zero 
and S(x) need not necessarily be a monotonically non- 


decreasing function . 

Ihe discrete size grinding parameters in 
Eqs. (5.27) and (5.28) become 

/ I'^(v) M(v 5 t) dv 


syt) = 


^ 1-1 


/ M (v, t ) dv 


X. 


(6.50) 
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selekage f imc'tion "blie foiaa is nscessai’ily 


I(y,x) = li(v) 


1*2 (2:) 


(6.27) 


wh. 62 ?G li(x) is ajiy posi'fci've definiiie func'tion over "tiie 
range 0 < x < , It follows that 



and from Eqs. (5.21) , (5.22) and (6.27) 
S(v) = h(v) 


and 


B(v,x) 


1*2 (^) 


(6.28) 


(6.29) 


It should be noted that in contrast to the RSP model , 
in this case S(0) need not necessarily be equal to zero 
and S(x) need not necessarily be a monotonically nor^- 
decreasing function . 

The discrete size grinding parameters in 
Eqs. (5.27) and (5.28) become 

\-1 

/ I'^('v) M(v, t) dv 

SjC-t) = (6.50) 

i-1 

/ M (v, t ) dv 
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Or, compactly 


where 


Sj_(t) = » ■fc) 

^i-1 

/ 1-1 (v) M(v,t) dT 


(6.51) 


P(x^t) 




^i-1 


(6.32) 


f M(v,t) dv 




Rate parameter is therefore , a function of time . ind 

lig (^i_l) - l2(x.) 


B. . 


i^~rxT7 


(6.33) 


which is independent of tone and shnilar to . in the 
RSP model in Eq. (6.22) . Substitutions in the'master 
equation Eq. (5,10) generate the size discretized SSE 
model for grinding kinetics 

dM. (t) 

= - I<2^^i) ^’(3:^,^) M^(t) 


,, i-1 


(6.34) 
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Incidsnijallyj “tliis is ilie first tijiie this model has been 
clearly delineated , She inner breakage parameter for the 
SSF model is obtained from Eq, (5.30) 


X 


1-1 


/ b-] (y) M(Y;t) dY 

tr ^ 


bi,i(t) - “ 


^i 




(6. 35) 


f S(v) M(v,t) dr;- 


^i 


The selekage function in Eq. (6,13), which is the 
product of the selection and breakage functions giYen in 
Eqs« (2,56) and (2.57), respectively,, constitutes a valid 
example of separable s^ekage function , In order to compute 
the time valuation of the rate parameter in the SSE model 
arising from the discretization step aM also the time 
variation of the inner breakage parameter, we substitute 
Eq. (2.59) into Eqs. (6,30) and (6,35) ? and obtain the follow... 
ing explicit expressions for S^(t) and b^ ^(t) 

a 


Ax. 


^i-1 


/ 

X, 


,.a-1 


[-jY ( yrry) 1*^ 


s^(t) 


a 


^i-1 p-1 

I T [-• si- ) I dv 

X^ 


( Sv30 


and 
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b. .(t) = 1 

X jX 


^i-1 a-1 a 

^ ^ C--a< ) ] ar 

Z' ^ r _1 ^ _v _ \“ -1 


, (6.37) 

lor simiilation, the yah/-,es of Zq , i. , and /i^iO) were 
k^t the same as in the previous section , Two combinations 
of a and p were chosen - (1) a = 0,5 and p = 2 and (ii) a = 1.5 
and p = 0,5 • In the first simulation experiment, computations 
were carried out for 6 = ‘^2 as well as 6 = ^2 , The particle 
size distribution cuives at the start of the experiment i.e. , 
at t = 0 and at the end i.e, , at t = 9 minutes, are shown in 
Figure 6,5 , Ihe results of the variation in S^(t) as a 
function of grinding time are shown in Figures 6.4 and 6,5 
and those of the time variation in b. .(t) are shown in 

X ^ x 

Figures 6,6 and 6,7 • Ihe rate parameter values have been 
normalized to S. (t)/s. (0) , It can be seen from these 

JL 

figures that the change is maximisn for size index i=1 and 
decreases progressively with increasing size index , This 
trend is quite expected because the particle size distribution 
in the coarse size range changes with time much faster than 
the size distribution in the fine size range , inother 
significant observation is that the time variation of 
discrete parameters is conparatively low when narrower sieve 
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6.3 Initial and final particle size distributions 
in simulation experiment for theSSFmodei 
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size intervals are used for describing the particle size 
distribution. 

It should be pointed out that , in practice, more 
pronounced variations maj be expected to occur because, 
unlike in this example, the feed size distribution cur/e 
in the coarse size range does not normally exhibit a 
pronounced curvature , Ibr illustration of this point, the 
values of the time invariant rate parameters (unifoamiLy 
distributed particle size distribution - Bass model } a:id 
those in the above simulation experiment at t = 9 minutes 
have been compared in Table 6.2 . The percent change in the 
values of rate parameters is now almost tv/o times of that 
shown in Mgures 6.4 and 6,5 . Many other simulation 
experiments with different combinations of a and p showed 
that if the difference in the values of a and p is increased, 
time variation of discrete size parameters becomes more 
significant , In one of such experiments with a = 0,5 and 
p = 4 , the per-cent increase in the values of the rate 
parameters was found to be approximately twice of that 
reported for a = 0,5 and p = 2 , 

The initial and final particle size distribution 
curves for a = 1,5 and p = 0,5 are shown in Figure 6,8 , 

The time variation of the rate and inner breakage parameters 
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TABLE ~ 6.2 

Comparison of the values of time invariant rate parameters 
anl the rate parameters at t = 9 min, in the simulation 
experiment (Eig. 6. 3, a = 0.5 , p = 2 ). 


1 


1 

3 

6 


6 = f 2 

\ 

6 = 



S. (time 
^ inv, ) 

S^.(t=9 
min. ) 

Percent 

difference 
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S. (time 
^ inv. ) 
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1 
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is shown in Figures 6.9 and 6.10 , respectiveLy . In contrast 
io "blie Figures 6.4 and 6,6 where rate parameters increase 
with grinding time, Figure 6.9 shows that the rate parameters 
constantly decrease with increasing grinding time . In fact, 
simulation experiments showred that this is generally true 
for all the cases where^is greater than p , Similarly , for 
a < p , the rate parameters were always found to increase 
with increasing grinding time . Ihis observation can be 
explained in the following manner . For the case a > p , 

(v) is a monotonically increasing function in the particle 
size V , As grinding proceeds, the relative mass density of 
particles [ M(v, t)/l!ilj^(t) ] decreases for the coarser particles 
anh it increases for the finer ones , It, therefore, follows 
from Eq. (6,30) that due to the increasing proportion of 
finer particles in the discrete size interval, whose 
contribution to the rate parameter (i. e. , l^(v)) is relatively 
small, the rate parameter decreases with increasing grmiding 
time . Similarly, it can be argued for the case a < p that 
the rate parameter always increases with increasing grinding 
time . It will be further seen, that the logarithmic - linear 
plots of M^(t) vs. t for a single size material always tend 
to concave downwards for a < p and they concave i:pwards for 
a > p . Examples of both the cases can be found in the 
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Fig, 6. 8 Initial and final particle size distributions 
in simulation experiment for the SSF model 
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grinding literature (41,82,84! . 

finally 5 in view of the fact that the time variation 
is more pronounced in case of the inner breakage parameters 
than that observed for the rate parameters, ifc may be 
suggested that of the two discrete size grinding models, the 
model 'A* is a better choice, if the time invariant parameter 
models are to be used . 



141 


OHAPIER 7 

iiTJiiiisis OF sam spegim asiimim pimmmA 

In liglit of the thj?ee specialized size discretized 
models available to us, in the present ch^ter an attenpt 
is made to seek unatibiguous answers to the problans concerning 
the occurence of various ^ecial grinding phenomena that have 
been observed ei2?)irically. In particular , we shall examine- 
(i) first order disappearance kinetics for a discrete size 
interval, (ii) zero order production of fines, (iii) Bosin- 
Ranmler kinetics, (iv) Alyavdin kinetics, (v) difference- 
similar breakage parameters and (vi) size variation of rate 
parameters according to a power law. 

7.1 Eirst Order lisappearance Kinetics for a liscrete 
Size Interval 

In the literature, there has been an implicit or 
explicit assuHption in- treatment of discrete size grinding 
models to the effect that if the material breaks in the 
manner of a first order kinetic law, then the rate of 
dis^pearance of solids from a discrete size inteirval will 
also conform to the same law I 72 j . ^.orrigtimes, it is 
stated that this is a ^pTOximate hypothesis bub, as Harris 
j 45 I has pointed put, the approximations which lead to these, 
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kinetics for discrete size intervals have never teen clearly 
defined. Examination of Eq. (5.26) shows that the expression 
for the pure dis^pearance kinetics is 

dMj_(t) 

- - f I-(v,Xj_) M(v,t) dv (7.1) 

1 

It can he seen from Eq. (7.1) that the rate of change of 
Mj_(t) cannot he explicitly expressed as a function of iL (t) 
only. Thusj in general, it is not possible to realize any 
order kinetics, in terms of l/^(t) , leave aside the first 
order law . However it will he seen that the restrictions on 
the particle mass density function M(x,t) in Eqs. (6.1) aid 
(6,2) , resulting in the Bass model, and/or the restriction 
on the selvage function in Eq. (6.15), resulting in the RSP 
model, are sufficient conditions for first order hypothesis 
to prevail. But it shoiild he noted that the rate parameteis 
in these two models, given in Eqs. (6,3) and (6,21), have 
couplet ely different forms and structures. 

It should he pointed out that the first condition of 
uniform size distribution, within individual size intervsls 
can he met by only extremely narrow size intervi^s, which 
cannot be satisfied in practice . SDhe second condition which 
though represents a highly unlik^y situation, provides 
iEpoitant clues to this pro hi It will he sem from 
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Eq. ( 5 . 27 ) tnat lesser is the dependence of function L on the 
parent particle size v, greater is the size interval width * 
that can be afforded in practice to obtain a reasonable 
first-order- kinetics in the conventional sense. Of coarse, 
when the selekage function is conpletely independent of the 
size of the parent particle, the width of the discrete size 
interval is immaterial. 

Austin I 28 j has pointed out that when a discrete 

size interval is too wide, the conventional first-order 

disappearance kinetics plots are either concave upwards or 

concave downwards* He has ascribed this phenomenon to a 

number of factors, which have been mentioned earlier in 

Gh^ter 2. A close infection of Eq, (5.27) for the rate 

parameter j however, reveals that the answer to this 

problem lies, in fact, with the Selekage function L(v,x). Eor 

a discrete size fraction i, l(v,Xj_) , which represents the 

fraction of material of size v which falls out of size 

interval i duiring the breakage, is normally e35pected to be 

either a monotonically increasing function in v or a 

monotonically decreasing function in v , over the size range 

X.* < V < X. * Moreover, as the grinding proceeds, the 
1 i— 1 

proportion of coarser particles in the discrete size interval 
decreases and that of the finer ones increases. Hence, it 
is clear from Eq. (5.27) that when l(v,x^) is a monotonically 
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increasing function in v, lesser proportion of material of 
size interval i vail break out of it at liigber grinding times. 
HenoSs Sj|_ will exhibit a decreasing trend with increasing 
grinding time, and the dis^pearance kinetics plot will 
therefore, concave upwards, The reverse holds when l(v, x^) 
is a 300 notonic ally decrrasing function in particle size v. 
Obviously, depending on the initial size distribution of the 
particles, and the nature of the function l(v,x^) over the 
size interval x^^ < v < x^_^ , even more conplex kinetics may 
be observed, 

finally, the other factors (Gh^ter 2,5) mentioned by 
iustin ( 28 I , may augnent or nullify the effects of 
initial particle size distribution and the selekage function 
on the dis^pearance kinetics , However, in view of the 
results obtained on time variation of the rate parameters 
in CShapter 6 and those reported later in Oh^ter 9, it would 
seem that these factors are of secondary iii 5 )ortance and 
the main deteimining factor is the seldcage function. 


7.2 Zero Order Production of Pines • 

We let X = Xj^ , and rewrite Eq., (2,34) ^ follows 


dP. (t) 
dt 


^0 

/ L(v,x.) 

^i. 


6v 


dv 


C7.2) 
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The phenomenon of the zero order production of fines requires 
that ri^t side of the equality in Eq. (7.2) be a constant, 
at least for some initial period of time. ¥e shall now 
scrutinise the discrete size models available to us for 
possible justification for this phenomenon. 

In case of the Sass model we first write Eq. (7,2) as 

i ^i 

"■"dt ^ Lf S(v) B(v,x) to] M(v,t) dv (7.3) 

and then invoke the stipulation of uniform distribution 
within the individual size intervals to obtain 


dP^(t) 


i M.(t) ^i 

^ 1 n ^(v»x) dv (to 

iJ""* DO X. 

^ ( 7 . 4 ) 


dP^.(t) 


i M.(t) ^ 

S 2 / 

j=1 k>i+1 


>-1 ^ j -1 

Z / / S(v)B(v,x)dv dx 


( 7 . 5 ) 


Substitution of Eq. (2.40) gives 


dE^(t) 

dt 


Z T ] 


k>i+1 


(7.6) 


Por the RSP model , we substitute Eq. (6.16) into 
(7.1) and recall that P( 2 q, t) = 1 for t > 0, hence 
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we have 
dP. (t) 

— r 1 “ ] ( 7 , 7 ) 

Por the SSP model, we sub^itirbe Eq. (6.26) into 
Eq, (7*1) aiid obtain the following expression 

dP.(t) ^0 

— 51 - M(v,t) dv (7,8) 

1 

which we write as 

ai'.C't) , i ^3-1 

- d^— = J' M(T,t) dT (7.9) 

3-1 Xj 

Substitution of Bj,. (6.30) followed by substitution of 
Eq. (6,35) gives 

dP (t) i 

Equation (7.10) is apparently same as the one presented by 
Herbst and Puerstenau ) 58 | for the rate of production of 
fines. However, it should be noted that in contrast to their 
ass uopt ions of n eligible inner breakage parameter and 
environment dependent rate parameter, we have derived Eq, (7,10) 
by discretization of an inherently linear gcinding mod^; 
the time dependent rate and time-invariant bre^age parameters 
are the characteristics of the SSP model. 
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Inspection of the three expressions in Eqs. (7,6), (7,';' 
and (7»10) for the rate of production of fines shows that the 
most natural candidate for the realization of the zero oid.er 
production phenomenon is the fiSf model. We ai’gue that if 
is much smaller than the smallest particles in the feed then, 
at least for initial period, (t) << 1, and 


dP^(t) 




(7.11) 


Ooii?)aring Eq,. (7.11) with Eq., (?.. 116) atid (2.124), it can he 
seen that the function is nothing hut the function I2 » 
Substitution of Eq. (6,21) into Eq. (7,11) gives 


dE^(t) 


^ S. 


( 7 , 12 ) 


dt i 

Substitution of Eq, (6,16) in Eq. (5,29) gives for the RSF 
model 






S. 

1 


(7.13) 


Erom which it follows that 




B. . 

is3 


(7.14) 


Combining Eq. (7.14) with Eq, (7.12) gives the restrictive 
interrelationship between the rate and breakage parameters 
proposed by Mika et al, j 42 | , and Herbst and Euerstenai l58{ 
as a sufficient condition for zero order production of fines 


dE^(t) 

di 


S. B. . 

3 


(7.15) 



1^',8 

This condition is then a conseq.uenoe of the SSP model. 

It should te pointed out that only when we specialize 
the function = A z“, the Eq. (7.12) becomes 

at - A (7.16) 

?/hich is the form observed by Arbiter and Bhrany | 83 J . 
Since ' there is no reason why should always be of the 

form A , for this reason, it is possible that in certain 
grinding systems the zero order production rate constants 
may not conform to the form x? . 

7,3 The Ro sin-Ramml er Grinding Kinetics 

Starting from Eq. (2,106) we shall derive an exact 
solution for the JiSE model. Substitutions of Eqs, (6.21) and 
(7,14) causes the summation term to vanish, and a straight 
forward integration of the remaining equation gives 

\(t) = R^(0) exp [ - l 2 (x.) t 2 . C7.17) 

This is the discrete size analog of the continuous size 

Ro sin-Ramml er type distribution in Eq. (2.54), When the 

# 

oc 

reduced selekage function is equal to A x , we have the 
Rosin Rammler distribution 

\(t) = E^(0) ejp[-Ax^“ t]' (7. IS) 
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We like to stress that the BSI model holds when a plot 
of In against time t is a straight line. On the other 

hand, a plot of In In [ / E^(t) ] against In x. will 

be a straight line (with slope equal to a ) only if IgCx) = 

Ax •* » 

Over the years, the conminutionists have developed 
an uncertain altitude towards the iiosin - Rammler distribution. 
Many workers feel that by virtue of its relatively sinple 
form, the Rosin-Rammler representation of the grinding kinetics 
cannot conceivably do justice to this coii{)lex system, lonethe- 
less, the Rosin-Rammler kinetics continues to attract a steady 
patronage, either overtly, as for example, by Austin and 
luckie j 105 1 , Ruruya et al, | 106 [, and Gupta and K^ur 
I 19 ( , or covertly, as for example, by Schonert | 36 j , 
and Herbst and Ruerstenau j 58 | , 

7,4 The Alyavdin Kinetics 

Harris aid Ghakrdvarti | 46 j have shown that the 
Alyavdin size distribution function fits a variety of 
grinding data quite well, ^ecially where the Rosin-Rammler 
distribution turns out to be inadequate . The ilyavdin 
distribution is 

a(x,t) = E(x,0) a® (.'f'(x) 1:^+® ) , s. 0 (7.19), 
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where ^ is a function of particle size, and the e^onent 
a IS nominally a constant, although it may exhibit a weak 
d^endence on x and/or t . Harrib 1 38 | found it convmient 
to introduce the exponent a where dq^arture from the Eosin- 
Eammler t^e kinetics occurs . His derivation essentially 
PBcjuiicss "tiiOT Lrh.6 S6l6kag6 func'tioB be wpi’fc'tein es 

I,2(-'^,x,t) = t® (7.20) 

We shall present si up le arguments for finding an q)proximate 
solution to the SSH model. Prom Eq^s. (5.29), (6.26) and (6.31) 


Sj(t) 

Substitution in 



Eq. 


= s,(t) 




(2.106) gives 


t) 


( 7 . 21 ) 


dH.(t) i-l 

■sTTETd- = - 


E(x._^^,t) -P(x ,t) 
FCx-';-t) 


] RjCt) 

(7.22) 


If the function P is only weakly dependent on x such that 
the following inequality holds 


i-1 

Z 

d=1 



, t) - P(x. , t) 

piifTTr^ — 


H (t) 


1 « 1 


(7.23) 


then to a first q)proximation we may ne^ect the summation 
term, and solution to Iq. (7.22) is given by 

t 

Ej^(t) = B^(0) exp / P(Xj^,t’) dt' ] ,(7.2+) 
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where the integral term is a weak function of z . Equation 
(7.24) le the general form of .the Alyavdin distribution in 
Eq. (7.19) . Simulation studies have shown that the S(x,t) 
values computed from the similarity solutioh in Eq.(2.59) 
for sparable seldcage function in Iq. (6.13) (with a dose 
to p) , can be plotted as per the jSLyavdln equation and, in 
most instances, the deviations are well with In the 
experimental errors in the data (Figure 7.1). Ihe point to 
be stressed is that the Alyavdin distribution need not arise 
from any inherent time - dependence of the selection and 
breakage functions, as in Eq. (2.63) or Eq. (7.20), but it 
could well be simply an approximate solution to the integro 
differential equation of grinding in cumulative fraction 
retained mode for u class of time- invariant selekage functions, 
as the one given in Eq. (6.15) . 

7.5 Difference - Similar Breakage Parameters 

Eor cumulative breakage parameter B. . , the difference 

t) 

similar characteristics in5)ly that 

®l,d = h-d <7-25) . 

Inspection of Bqs. (5.28) and (5.29) indicates, however, 
that, in general, the breakage parameters vary with time . 
Eurttiermore, even if the cumulative breakage function is 
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normalizable in following manner 

B(v,x) = B(x/v) (7.26) 

in the general case we are unable to attribute no rmaliz ability 
i, e. tile difference similar form to tbe cumulative breaica,ge 
p aramet er, 

Brom Eq, (5.29) tbe expressions for tbe cumulative 
breakage parameter in three discrete size models delineated 
previously are as follows. Bor tbe Bass model 

X. 




I 

X . 

■— f] 




l(v,x^) dv 


/ Ii(v,x.) dv 


(7.27) 


and for tbe RSB and SSB models 

igCxj) 




(7.28) 


In case of tbe Bass model, if and only if, tbe selection and 
breakage functions are given by 

S(v) = A v“ (2.56) 

B(v,x) = B(x/v) (7.26) 

tbe breakage parameter in Bq, (7.27) is differoace- similar 
(noimializable) as follows 
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if and only if, the selection and breakage functions are 
given by ©js* (2,56) and (7* 26), tbe rate parameter is 
given by following power law 

A X ” ■} 

% ^(1,6 ) (7.33) 

Similarly, in case of BSP model, if and only if, the 

oc 

selekage function is given by l(v,x) = A x , the rate 
parameter obeys the pow/er law and the expression from Eq. 
(6.21) is 

a 

Si = Ax^ (7.54) 

It should be pointed out that except for SSP model, 
the phenomenon of difference - similar breakage parameter 
is always coupled with the phenomenon of a power law for size 
j variation of the rate parameter. 

i’ 

‘i 

I 



CaapiER 8 


IDMIIPIOillON OP IHE GRINDIIG SISIM 
IF TERMS OF A Mlli FUFCffilOF 


As pointed out earlier in C3aapter 2 , there does not 
exist at present a practical and reliable method for estimation 
of the two basic functions S(x) and B(v,x) , which uniquely 
characterize the mill-material grinding system , ilthou^ , 
a number of approximate methods have been developed for the 
estimation of the size discretized paraneters y a unique 
identification of the grinding system is still not possible 
because of the d^endence of these parameters on the particle 
size spectra . It may then be concluded that at presoat no 
convenient scheme is available for a true identification of 
the grinding systan • In this chapter , a novel approach is 
presented for identification in terms of a mill function , 
which is based on a continuous size and discrete time 
representation of the grinding systea , In the next chapter 
this approach has been further extended for estimation of the 
size discretized rate and breakage paraaeters in the model 
»A‘ in Eq. ( 2.79 ) . 



157 


8.1 Poimulation of the G-rinding Model 

It hss hesH i?6cogiiis6d that the 3 ?e exist thD?ee 
mechahisnis for* size a?edu.ctioii of particles , nsmely , impact 
fracture , chipping and abrasion { iCf7 [ , ilhere are unfortuna- 
tely conceptual difficulties in arriving at unanbigLous 
definitions of these mechanisms . We prefer to define impact 
fracture as catastrophic failure of particles resulting in a 
number of fragments of substantially smaller sizes than the 
parent particle . Chipping entails detachment of small fragments 
from corners and edges , the resulting size i^ectrum of the 
progeny is discontinuous in the sense that it comprises one 
particle of size comparable to the parent particle and a ntmber 
of fine fragments • Por a sharp distinction between the 
chipping and abrasion mechanisms it is desirable to concq>tual- 
ize that in the latter phenomenon reduction in particle size 
occurs by continuous removatl of fine matter from its surface 
due to rubbing and friction , In theory at least , if impact 
fracture and chipping are the events operating and only a 
fraction of material is selected for breakage by these 
mechanisms at any instant of time , then no matter how long 
the material is ground t some amount will always r^ain 
unground • On the other hand when abrasion mechanism also 
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prevails then, by definition , all the particles in the 
initial feed must undergo seme size reduction even when the 
grinding time is very short . Neither the continuous variable 
grinding model in E^. (2.36) and , therefore , nor its discrete 
size versions in Eqs. (2.79) and (2.86 ) take cognizance of 
the abrasion mechanism which should contribute not insignific- 
antly to the overall size reduction process when relatively 
soft materials are ground , e.g. limestone . 

In our approach we consider a grinding system in 
which the particulate charge undergoes size reduction by 
fracture , chipping , as well as , abrasion . It is assumed 
that the particle population is sufficiently large such that 
meaningful statistical averages may be defined for the various 
functions onployed in formulation of the grinding model . 
Eurthermore , we shall be concerned with a system which is 
time-invariant and is also linear in the particle mass 
fractions M(x, t) . We next define a mill function in the 
frequency form as e(v,x,t), x < v , where 6(v,x,t) dz is the 
fraction of ground particles in the size range x to x + dx 
when an unit impulse of size v is ground- for a fixed time 
interval t , dearly 6 is a nornalized frequency function , 
hence , 

V = 

/ ©(v,x,t) dx *= 1 

0 


( 8 . 1 ) 
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Tile mill function in form of a r>nmnio+-?^ x •, x. 

xuiiu ai a cumulative distribution function 

can be written as 


e (v,x,t) - / e(v,w,t) dw (8.2) 

Por a set of batch grinding experiments equally 
spaced in time , where each time interval is of duration t, 
let the successive time integrals be denoted by k=1, 2, 3, . 
and the grinding time at the e n d of k-th time interval by 
tjj. « Applying mass balance on the particles of each size , we 
can directly write^the following expressions for the output 
of the mill after t grinding time ' 


^0 


M^(t^) - - / M(v,t^^^) e(v,z^ , t) dv , i=1 , 


k=1, 2, 3, ... (8.3) 


and for i > 1 , we have 


^i-1 


:^(tj^) = 9(v,x. , t) dv 


+ 


/ M(v,tj^_^) [ 0(v,x^__^ , t) - e(v,x^ ,t) ] dv 

^i-1 

( 8 . 4 ) 
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Substitirbion of Eq. (2.55a) in Eqs. (8.3) and (8.4) , followed 
by integration by parts gives 


M^(tj^) = 1 - e(xQ , 


■fc) + / e(v,x^,t) dv 


(8.5) 


and for i > 1 


- 1 - + 9(x^ , , t) _ , X., t) 

^0 

+ / ^’(v,t^_^) e(v,x^ , t) dv 

i 

^0 

- / , t) dv (8.6) 

^i-1 

Since M(v,tj^) is not directly accessible frcm size distribution 
data , the advantage of Eqs. (8.5) and (8.6) over Eqs. (8,3) 
and (8.4) lies in the fact that it is convenient to fit 
suitable orthogonal polyncanials to E(x^ , t^) in x . The 
resulting expression is then substituted in Eqs, (8.5) and 
(8,6) for an a-ialytical evaluation of the integral terms , 

It is obvious that the knowledge of the function Q 
is sufficient for identification of the systan model , which 
is now discrete in time . Ihe function 9 c is subject to a 
number of constraints t 
(i) 0(v, 0, t) — 0 , for all v 



161 


(ii) e(v,-v,t) = 1 , for all v 

(iii) 0(v,Xjt) is a moncftonically non-decreasing function 
in X , i. e. , e(v,x, ;t) > 0 , for x < v . 

8.2 DetenQination of the Mill Function and Dsaon strati on of 
the Proposed Identification Schone 

For a danonstrat ion of the proposed method of 
identification , we have chosen Berlioz grinding data | 108 j 
for dolomite ground dry in a hatch hall mill with particle 
loads of 3500 gm. and 2640 gm. . In absence of any ‘a priori* 
knowledge of the functional foim of 0 we assimed the. following 
generalized polynomial expression 

_ = 2 2 3 

e(v,x,t) = aQ+a^x+a 2 V+a^x +a^ 3 cw+a^T +agX 

+ ay + agXY^+ ... (8.7) 

It chn easily he shown that under the first two constraints 
mentioned ahoye Eq. (8.7) reduces to 

0(Tr,x,t) = ( I ) + xCv-x) [ bg + + bgv + b^i 

+ b^xv + b^T^ + 1 (8.8) 

where 

bo =- 


bi % 
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bg = - a? 

........ and so on . 

Wi-thoat violating the first two constraints , Bq.(8.8) 
can be modified to perhaps a more versatile form 

= ^ Pi P2 P3 P4 P4 

e(vjX,t) = (- ) +x V (v - X ) 

X [■bQ+b^x+b2V+b^x^+b^xv+b^v^+ ... ] (8.9) 

where the exponents ? P2 P4 greater than zero . It 
is noted that Eg, (8.8) is a i^ecial case of Eg. (8,9) when 
Pl “ P2 ~ ^4 ~ ^ ■■ P3 ~ ® coefficients b*s and 

P-j ? P2 * ^3 ^4 determined from experimental data 

minimizing a least sguares type error function 

D ^ 

3 * 1 2 ^ 2 

Er = I E (8.10) 

lata for three grinding time inteirvals t^ = 60 , 80 , and 100 
revolutions was employed (t — 20 mill rev.) . A combination 
of random search method | 109— 115 1 initial search 
followed by Rosenbrock method ] 11 6-1 .19 I for fine convergence 
was used for the non-linear optimization . Initial trials 
showed thst in the polynomial in Eg, (8,9) the constant b^ , 
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■fclie coefficients of pure x 'terms , the cross product terms 5 
and the terms of order' more than 5 were about 3 orders of 
magnitude smaller than the coefficients of first five pure 
V teims • Elimination of these terms in the polynomial did 
not increase the value of the error function Er by more than 
2 pe2>-cent . Consequently , the form of the function 0 was 

modified to 

2 - h ^2 P 3 ^4 ^4 

0(v,x,t) = (“) +3: V (v - x) 

2^4-5 

X [c^v+ C 2 V + C^V^H- c^v + c^r ] (8.11) 


which had only 9 paraneters to be estimated . !Dhe values of 
these parameters , aid the experimental ard computed size 
distributions for 3300 gn, load of particulate charge are 
shown in Tables 1 and 2 , respectively . It can be seen that 
the predictions of the size distributions are very good , 
Indeed j in comparison with the large number of simulation 
stt Qrpt s reported in the literature j this sesns to be the 
closest that has ever been achieved . 

In order to reduce the number of parameters in 
function 8 , a 7 parameter form was next tried 


Pl ^2 

0(v,x,t) “ ( I ) + C ( I ) 



3 


X c d^v + dgT + 



( 8 . 12 ) 
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3)i£LE - 8. 1 


Estimated values of the parameters in Eq.(8.11) for 5300 gn 
of parfcicjle load. 


Parameter 

Value 


43.827177 

P2 

00.890110 


-43.721428 

h 

42.835910 


0.482328 

02 

-1.660828 


2.179485 


-0.373028 


-0. 378885 
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Ii£LE - 8, a 

Gompa.Z'ison of Experimenl/al and Oompufed Size DLsfribu'fcions, 
Berlioz Bata for 3500 gm. Load of Bolcmite 1108} 
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Tlie estimated values of the parameters for 3300 ga. as well 
as 2640 gm, load of particulate charge are given in lable 3 
The predicted size distributions for 3300 gn. load were now 
within + 1 peiv-cent of those reported in Table 2 , The 
experimental and computed size distributions for 26-1-0 gm. 
load are shown in Table 4 * -^gain, the agreement is quite 
good . 

Por both the sets of data j four to five widely 
different starting guesses for the parameters in the three 
polynomials in Eqs. (8.10) , (8.11) and (8.12) were tried . 

The final parameter estimates did not show more than + 1 per- 
cent variation , The variation in the value of error function 
was also within + 2 per-cent , Port her evidences in support 
of the validity and the accuracy of the est'imated mill function 
shall be given in the next ch^ter . 

Pigures 8,1 and 8,2 show the nature of 0(v,x,t) 

* curves plotted as a function of dimensionless size x/x^ on 
abscissa scale , for feed impulses of size v located at 
various mesh openings v = x^ . In most favourable case when 
the impulse is located at 10 mesh opening (v=1 680 microns) 
about 18 per-cent material only is ground finer than the 
adjacent 14 mesh (x=1190 microns) in 20 mill revolutions . 
Indeed, more than 75 per-cent product' lies in the size range 
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lABiE - 8.3 


B^tlnatea value, of the pa.»ete.e In (e. 12 ) f„, 
3500 gt. and 2640 gm. partioie loads . 


Parameter 

__.Humefioal Value 

3300 021 . particle 

1 oad 

2640 gm. particle 
load 

Pi 

43.732543 

37. 279627 

P2 

0. 907984 

0,827502 

^1 

0.482858 

0.469503 

d 2 

-1.660334 

-1.587290 


2.179932 

2,202055 


~0. 372955 

-0.385656 

% 

-0.378812 

-0.412245 
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lAELE - 8.4 

Comparison of Experimental and Computed Size Distributions, 
Berlioz Data for 2640 ga. load of Dolomite ] i08 t 



Size 

index 


60.05 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 
12 


13.59 

8.54 

4.87 

3.13 

2.23 

1.82 

1.40 

0.94 

0.74 

0.56 

0.54 


60 rev. 


80 rev. 


100 rev 


Exptl. 

Predicted Ixptl. 

46.88 

46.90 

36,56 

15.43 

15.40 

15.82 

11.01 

10.68 

12,71 

6.73 

6.94 

8.21 

4.66 

4.75 . 

6. 06 

3.38 

3.48 

4.42 

2.84 

2.78 

3.84 

2.09 

2.15 

2.81 

1.40 

1.52 

1.87 

1.09 

1.19 

1.48 

0.79 

0.90 

1.10 

0.94 

0.86 

1.14 


Predicted 

E3q)tl. 

Predicted 

36. 63 

28.11 

28.61 

15.80 

15.41 

15.33 

12.26 

13.65 

13.35 

8.69 

9 . 42 ' 

10.13 

6.18 

7.35 

7.41 

4.60 

5.47 

5.59 

3.65 ^ 

4.85 

4.44 

2.84 

3.55 

3.48 

2.08 

2.26 

2.60 

1.62 

2.01 

2.02 

1.24 

1,42 

1.57 

1.11 

1.52 

1.36 





X/Xj 


Fig •8.2 Cumulative mill function as a function of scaled size 
x/xi for feed impulses of size v located at various mesh 
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to 0.9 . This strongly indicates that in course of the 

size reduction , chipping and abrasion played a significant 
role . The effect is even more pronounced when the feed 
impulse is located at smaller sizes . 

It is interesting to note that in our choice of the 
function 6 in Eqs. (8.11) and (8,12) we could not guarantee 
the third constraint regarding the monotonic nondecreasing 
nature of this function , However , the computed curves show' 
that this condition is not violated , 

Figures (8,3) and (8,4) show 9(v,x^ , t) curves , 
that is, the fraction of feed material ground out of a sieve 
size interval , when the feed impulse ranges over two 
adjacent mesh openings , ^ all cases , w'hen 

the impulse is located anyi/vhere in the size range £’ + 0.2 
^Tc. ^ - X.) 1 to X. ^ , most of the material after 20 mill 
revolutions is still coarser in size than the 10 v?er mesh 
opening x^^ . In other words , a large proportion of solids 
after size reduction still remains within the original size 
interval . This result confirms that the inner breakage 
parameter b- . is not negligible , as was pointed out 
earlier in chapter 6 also , 



330Cgm particle load 



• ( i ^ / V ^ j, ^ P 

Cj'Tiuiative mili function for x=xj as a function of(v-xi)/(Xj„. 5 -xj)^ 
representing the fraction of feed material ground out of sieve 
size interval t when the feed impulse ranges over the Size 
ranne xi<v <'x. i 33n0qm particle load- 
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8.5 Conclusions 

A batch ball mill gciuding system has been identified 
in teims of a mill function -Q (T,Xjt) using a continuous size 
and discrete time mathematical model , Ihe procedure presented 
here sears to be of reasonable generality , Ihe role of 
chipping and abrasion mechanisms has also been properly 
delineated . 

The proposed identification scheme has the follov/ing 
advantages s- 

(i) The mill function provides more insight into the 
breakage behaviour of the particles than that can be 
obtained from the size discretized paraneters . 

(ii) In contrast to the size discretized parameters , 
which do not permit a unique identification of the 
grinding system because of their dependaice on the 
particle spbctra , mill functibh provides a more 
sound ard uniform basis for conparison of different 
grinding systems , 

(iii) The task of parameter estimation is considerably 
simplified when the mill function la used for identlfl 
cation. Eofe.exa.ple , In case of the two grinding 
systaEs considered above , only 7-9 unknown constants 
were retjuired to be estimated • Whereas . for 
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discrete size intervals, 78 time dependent size 
discretized parameters must be estimated 

(iv) Ihe proposed identification scheme is economical from 
the point of view of experimental labour , too« Only 
two to four batch grinding experiments . equally 
spaced in time , are sufficient for the identification 
of the grinding system , 

(v) The discrete time mill function can be anpolyed for 
estimation of S- and B. parameters in the more 

i j-j J 

commonly used continuous time and discrete size 
grinding model 'A' in Eq, (2., 79) ? 
next chapter . 


as shown in the 



OHiJTER 9 


ESTBL 


OP ME BaEiZj.GE P.JliI.;iE!rERS 


In absetice of any tlieoretical basis for prediction of 
n (n + 1) / 2 time varying rate and breakage parameters in 
Eq_. C5»10) , the paraneter estimation from noise corrupted 
experimental data for n as large as 10 or more, is not ati 
easy task . During the course of present investigation ^ a 
number of methods based on quasi -linearization j 120 , 121 ), 
Kalman filter ]122 , 123 j and also parameter optimization 
using polynomial approximation for time variation wore tried, 
which resulted in little success . It is then not surprising 
that all previous work in the area of parameter ostimation 
has been based on simplified models of the grinding kinetics, 
as reviowod earlier in Ghspter 2. In this chepter, our 
objective is to estimate these parameters using minimum 
possible simplifications . Since, we have had little success 
in tracking the grinding parameters in time, i^ order to 
malcc the problem tractable we have assumed that for 
appacopriate batch grinding data obtained over short grinding 
time intervals of duration t , the grinding parameters are 
sensibly constant in any given time interval k in the tine 
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range ^ \ • Hence the grinding nod el beccnes 

-grp- - - S^(k} M. (t) + s S.(k) B. .(k) n,(t) , 

i'~*l 52 j iij 

k-^1s2^ ••• ^ (9#1) 

Vi i i \ 

The cstination procedure presented here is based on 
mill function 0 ^ which has been defined end estineted 
in the previous Oheptcr , 

9,1 Devclopnont of the Proposed Scheme 

Prom Bq. (2.103) it follows that solution to 
Eq. (9.1) is 

M(tk} = 0(k) (9.2) 

where M(tj^) is the n x 1 vector of the product size- 
distribution, M(tj^__^) is the feed size distribution 
vector at the start of the k-th grinding intervsil , 

0 (k) is an xn lower triangular matrix operator which 
n^s the feed vector to product vector for k-th tine 
interval , The elenent o. .(k) of the natris: 0 is the 
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fraction of feed in j-tli size interval that is f c und in 
i-th size interval in the product after grinding for t 
time duration . It can be seen that the o-th vector 


in 0 matrix can be directly obtained by grindinv a 


single size feed of size j for a time period of duration 
t . jm alternative method of calculating 0 matrix, vdiich 
is based on the mill function, is given belov/ . 

By definition of e(v,x,t), we have 


= 1 


/ 9(v,x^,t) dv 

fi 

^i-1 

/ dv 

Xf 


(9.3) 


and 


o (ic) = 

X j J 


^ - = - = 

^ M(v,tj^_^) [e(v,x^_^,t)-9(v,x^,t) ] dv 


""j-l 


/ M(v,t^_^) dv 




(9.4) 


Substituting for M(v,tj^___^ ) from Eq. (2.35a), Eqs.(9,3) 
and (9.4) bocone 


^ L ® ^^i-1 » \\ 3 


X 


1-1 


+ / e (v,x^ , t) 


dv 


^i 


0 



Vi 


TVTT^.'IT 


(9.5) 
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X, j L® (■’'j-1’ ^i_i > ) 

-S(2^._^, , t) ] 


- I ^i-r “■ 1 

""d-i 

- S ^'('^»\>l) [ 6(v» » ■•= ) 

= 1 

-e(v, x^ , t) ] dv ; (9.6) 


Substituting the expressions for B and 5'j|_(k) in Bqs.(9.5) 
and (9.6) , 0 (k) matrices can be evaluated for each time 


interval k . 

Now, 

obtained for 


‘•V ‘ 

if the explicit alDeabric expressions can be 
each element o. .(k) in terms of the parameters 

J 


S^(k) and ^(k) , it shouldbe possible to back- calculate 

the parameters one by one . Neither Reid’s solution in 
Eq. (2.95) nor Herbst and Mika’s solution in Bq. (2.102) 
is convenient for this purpose . Ibis difficulty has been 
overcome by devising a ‘0,1’ variable suppression technique. 

Ihe procedure is as follows . 

It can easily be seen that the dia^nal elements 


\. 



°i,i simply exp (- t ) . Hence, all the n rate 
parameters can be first obtained immediately , 

1 close analysis of the structure of the 0 matrix 
reveals that if in a particular row one moves tovi/ards 
left, only one new breakage parameter is encounterc4 in 
every column . Moreover, the basic structure of the element 
o^ is as follows 








(9.7) 


It will be further seen that if all the parameters are 
known, substitution of 1 for and 0 for the rest of the 
M's in the Reid's solution to Mj;(t) in Eq. (2.95) would 
inmediatoly give us the value which will be equal to o. ^ . 

-L 9 J 

Since, we do iX't know B. • , we first calculate o. • as 

1 JF J ? J 

described above assigning B. . a value 0 . This gives us 

J-y J 

the constant m,. • in Eq. (9.7) . Repeating the same 

X y J 

calculation again, and now assigning B^^^. a value 1 , we 
get a numerical quantity which is equal to ( • 

Hence, using the known actual value of obtained from 

the mill function, and Eq. (9.7) > “t-be parameter 
be calculated . It follows that 
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In this way all the breakage parameters can be deteminea 
S6q.uent daily. It should be pointed out that unlike 
Austin end Luckie's approximate back calculation method 
j 78 I 5 back calculation procedure described above 
is exact and, therefore, no formal verification is as such 
necessary . It may, however, be mentioned that when tested 
against synthetic data, the above method did provide the 
exact estimates of the parameters , as expected , Por 
clear illustration of the procedure, a flow chart has been 
given in 4)pendix 3 . 

9.2 Estimation of lime Varying Parameters for 2wo Real 
Grinding Systons 

Since the estimates of the mill functions 9 and 

the orthog-onal p'jlynomials for Pj.(k) are already known 

for two of the batch grinding systems studied by Berlioz 

I -jos j > these may be substituted in Eqs. (9.5) and (9.6) 

in order to evaluate 0(k) matrices for grinding tine 

intervals 20—40 , 40-60, 60-80, 80-100 and 100-120 mill 

revolutions • Eor each time interval the pcirameters ^j_(b) 

and B. .(k) can be obtained using the exact back 
i > 3 

calculation method described above. 
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Figures 9.1 and 9,2 show the variations in the 
rate parameters with time for 3300 gn. and 2640 gn. loads 
of dolomite . Figures 9.3 and 9.4 show plots of some of 
the breakage parameters as a function of time • We note 
that manw” parameters ^ especially the parameters pertaining 
to coarse size interval s> undergo a distinct change from 
first to fifth time interval . In general, the time 
dependence is more pronounced for- the rate parameters than 
for the breakage parameters . Ihe breakage parameters 
are not difference similar as widely assumed hitherto 
) 42, 58, 61, 64 , 72 j , Ihe breakage parameter B. . . 

J-“r * j 

is considerably on the higher side when tne size indez 
i is greater than 5 or so, This again reflects the 
predominance of chipping and abrasion mechanisms in the 
size re<luction of fine particles . 

9.3 Definition and Calculation of the Reference Parameters 

Since model parameters are function of particle 
size distribution within the discrete size intervals, it 
is essential to fix a standard reference for comparison 
with different grinding systans . Ihe most convenient 
choice is a feed of uniformly distributed size distribution, 









ikage parameter 



Grinding time tj^,mill/ev. 


Fig. 9-3 -Variation o1 breakage parameters over dis«jet 
' ' intervals of time for 3300 gm- particle load. 






■ ^ Grindincj tinic rev. 

Fig' 94 Variation of :)rc'-:k.-:ge parameters over discrete 
intervals of, time for 2640 gm- particle loaf. , , 
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ground for a fixed period of tiae t . por a imifornly 
distributed feed Eqs. (9.3) and (9.4) reduce to 


O . . = 1 

X,1 


1-1 ^ = 

/ 6(v,x. , t) dy 

X. ^ 

1 

^i-1 - ^i 


(9.9) 




,^D-1 - = _ = 

/ [ e(v,x^__^, t) - e(v,x^, t) ] dv 

X • 

D 




( 9 . 10 ) 


Thus , a;pplying the back calculation method described above 
to the 0 raatrix whose elements are given in Eqs. (9.9) 
and (9.10) , both sets of the reference parameters can 
bo obtrunod . 

The lower triangular parameter matrices in Tables 
9,1 nJid 9.2 give the values of the reference p-rameters 
for tW‘j of the- Berlioz's grinding systems. The diagonal 
elcnontc arc the rate parameters for respective size 
intervals i. Uhe off - diagonal elements in i-th mw ana 

j-th column are breakage parameters . 

Mika ct aL. | 42 j have r^orted the values of the 
feed size rate paraiaetors , estimated from log-linear 
disqopearaaoo kinetics plots, for these two grinding 



!> 

o 

u 


H 

H 

•H 



a 

CO 





<D 

-P 



<P 


II 

§ 


11 -p 

Pt 


Cj 


« r <. 

Ph 



0 


s 

C3U0 

cd 


H 




CP 


<D 

Pt 



P0[ 





H 

<D 

P 


H 

cd 


O 



P 

CO 





P 

Pi 


o 

0 



B 


ncJ 

0 

H 

cd 

O 

d 

• 

0^ 

P 

rH 

o 

CtJ 

1 

% 

PI 

O 


ttO 

P«q 

a 

cd 

o 

•H 


CO 

P 


CO 

» 

EH 

' 

Cm 


B 

O 


CP 



4J) 



CO 

PI 


(1 

cd 



to 


OjO 

Pt 


PI 

0 


•H 

4^ 


'*0 

0 


Pi 

<^.-1 

s 


U 

Sh 


o 

cd 



o, 




Reference Parameters for Berlioz Grinding System, 2640 gai. Load of Dolomite \ 108 |, t = 20 mill 
Diagonal Elements are rate parameters and off-diagonal Elements are Bredcage Perameters . 



0.009 6 0.0153 0.0139 0.0145 0.0175 0 . j22 0 0.0269 0.0315 0.0368 0.0091 0,9037 0,2124 
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systems. Ttiey have also reported the values of the 
araiueters B. ^ , which were determined from estimated 
values of constants and Eci. (2.109) . 1 comparison 
uf these parameter estimates with the reference parameters 
determined above, is presented in Table 9.3 . It can 
be seen that except for aid , in both cases 

the agrccuait is quite satisfactory. It should be noted 
that in case of first two breakage parameters , Mika et al. 

T-n lic’ted that there existed inherent uncertainties 
in tliu procedure adopted by theni. Moreover , it should be 
pointed out tViot Berlioz's data does not oonEorm to 
aosta - type ktaetlos and henoe, strictly speaking, 

Eq. (2.109) ic not applicable . 
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TABIE 9.3 


Conpexisan of the parameter estimates rq,ortea by 
Miia et al. | 42 | and the reference paramet. 

obtained in the present study. 


'ers 


Paranet crs 

3300 gii. Pc 
lo 

.rticle 

ad 

2640 gm. Particle 
load 

i*ika et al' 
Est imat es 

s Ref. 
Parade 
ters 

1 Mika et al 
i Estimates 

j 

's Reference 
Parameters 

[ j.iin."' *J 

‘ 0.6210 

0. 61 28 

j 0. 6858 

0. 6673 

■^2,1 

0.406 

0.473 

0.389 

0.473 

^3,1 

0. 205 

0.142 

0.206 

0.130 

^4,1 

0.098 

0.106 

0. 1 09 

0.102 

^5,1 

0.063 

0.076 

0.062 

0.075 

\l 

0.044 

0.055 

0.048 

0.056 

^,1 

0.035 

0.040 

0.035 

0.041 

®8,1 

0.023 

0.029 

0,027 

0.031 

Bq . 
y f * 

0.019 

0. 021 

0.017 

0.023 

^10,1 

0.013 

0,015 

0,016 

0.017 

B11,1 

0.011 

0,011 

0,012 

0.013 

^12,1 

t 

1 

0.010 

0.008 

i 

1 

0.009 

0.009 


GHftPTER 10 


HVIPIRlCiiL GORilELAiriOES POR ViiRIiiSIGR OR IHE RAIE PiiRiiluETERS 
WlSii PARTI nUjjiiTE M-hSS -Alip BaI.I Sli^Sl 

in view of tiie fdoi uRtit voiy licoie qu.antitative 
info unit ion is available at present on the breakage 
mechaniaiis operating in a ball mill , it would be useful , 
from the point of view of the simulation studies , if 
sjstanatic procedures could be established for anpirical 
det eimination of some of the more important correlations for 
variation of the model parameters with the operating 
variables , In tliis chapter two such methods for deriving 
the appropriate correlations for the effects of particulate 
feed charge mass and the ball size on the rate parameters 
have been presented . It has been assimed that the Bass model 
is valid for the grinding systvsm 

lOol Effect of Particulate Mass 

As mentioned in Chapter 2 ^ it is nor 'a. priori' 
evident whether the influence of the particulate mass on 
the rate parameters is independent of particle size ^ In 
order to investigate this aspect we proceed by assuming that 
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it is so and the effect of the particulate mass can be 
incorporated in the following manner t 

S^CW) = H(¥) (10.1) 

where is the rate parameter for i-th size interval, 

is a characteristic constant associated with i-th size 
inteirval for a given mill and material , and H(¥) is a size 
invariant function of the particulate mass W . Our objective 
is to determine the function H(W) from experimental data , 
provided , of course , Bq. (10,1) is doaonstrably correct , 
We rewrite the discretized batch grinding equation 
in the cumulative retained mode R (2q, 2.106) as 

dR. (t) i-1 

CVl(W) 

- S.(W) B, . 1 iL(t) , 1 = 2, 3, ... (10.2) 

j -* J 

In Eq. (10.2) it has been assianed that B is independent of 
W |42 , 87 I • Combining Eqs. (10.1) and (10.2) we obtain 


a(w)'"k = [ - % ®i,3+1 ■ ^1,3^ ^ 


(10.3) 


We next define a scaled time 9^ as 


r = H(W) t 


(10.4) 
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Bq. (10.3) beccmes 


dRj_(:r ) 

“d 


i-1 

= - s. fi. (.r) + s 

^ ^ 3=1 



B. . . 

1 , 0+1 




(10.5) 


Tile solution to this equation can be fomally (vritten as 

Rj_( r ) = X . ( r ) (10.5a) 

It will be seen that if the time variable t is scaled by a 
factor H(W) , the basic foun of the equation governing the 
batch grinding kinetics as well as its solution in the 
domain of the scaled time variable 3" is invariant of W, It 
then follovi/s that the shape of tne ground particle size dis- 
tributions should be identical for a given degree of fineness, 
independent of the particulate mass W. In other words , the 
trajectory of particle i^ectra is invariant , only the 
grinding time t required to attain a specific distribution 
is scaled by the function H(W) . Berlioz data j 1081 for 
dolomite in Pigure 10.1 shows this clearly , even though 
the particulate mass and the grinding time differed 
considerably . Rather than depend on a chance of obtaining ; 
identical distributions from experimental data , we shall 
now develop a systoratic procedure for confirming tne 
validity of Eq. (10,1) , which concurrently deteimines the : 
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200 100 48 , 28 . . 1.4 

, particle’ SIZE , mesh 

Cumulative fraction finer than the statca S'ze 
for three values of the particulate mass demons 
trot ing the invariance . of the size spectra | 




nature of the function H(¥) . 

let r be the cimulatiTe fraction retained on -the i-th 
size after grinding time t^( r ) . From Eqs. (10.4) and (10,5a) 

H(W) t. (r) = H° t? (r ) (10.6) 

where the superscript denotes any arbitrary reference parti- 
culate mass , let 

(W) = H(W) / H° 



Thus , the ratio of grinding times to reach a giren cumulatiTe 
fraction retained for size index i gives the numeric 
value of H’ (Yif) . ii?om Eqs. (10,1) and (10.7) we obtain 

S.(W) = Hf(W) S? (10.8) 

With the knowledge of the function H'(W) and the rate 
parameters S? for the reference paiticulate mass , the values 
of the rate parameters can be obtained for any desired 
particulate mass W. 

Berlioz data jlOSj for 1320, 1980, 2640, 3300, 

3960, and 5280 gn, dolomite was plotted in the cumulative 
retained form as a function of grinding time . The reference 
mass chosen was 2640 ga, which resulted in near saturated 
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regime , Por eacii feed charge 30-35 ratios of t°(r)/t(r) v;ere 
computed from the graphs using all size intervals (i = 1-12) 
and different values of r in the range 0.3 to 0.95 . fflie mean 
and + standard deviation bound on H' (¥} are shovi/n iii 
Pigure 10,2 , It will be seen thot within the experimental 

and computational errors , - the latter due to graphical 
interpolation of data - , the values of H' (Yif) are reasonably 
constant in both the starved and over loaded regimes . 

Berlioz data for 660 gns. has been omitted frcm Pigure 10.2 
as it showed wide discrepancies in the values of H*(W) . 
Obviously , for a highly starved mill our basic assumptions 
in Eps. (10,1) and (10,2) are not valid , 

In order to find the variation of H*(¥), its inverse * 
is plotted as a function of particulate mass W in Pigure 10,2 . 
Por all practical purposes this curve can be approximated by a 
straight line , and the equation for H' (W) becomes 


where and k 2 are constants for a given material and mill. 
Combining Eqs. (10,8) and (10.9) gives 


Si(w) 




k-, + kg W 


(1C. 10) 


which is the required correlation for the variation of the 
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rate parameters wita the mass of the particulate material . 

In order to confirm the validity of Eq^. (10,9) , five 

rejDrer-ontotive sizes is plotted as a function of the scaled 
time T , \aiich was calculated from Eq.s. (10,4) and (10.9) . 
Since the scaling factor H° is a constant , it was arbitrarily 
set equal to unity for convenience . Eigure 10.3 shows that 
the restating curves are independent of the particulate mass, 

in agree . ciit Vvith Eq_. (10, 5a) . 

In Eq. (10.10) the value of the constant may turn 

out to uc of the same order as the tem kg W , For Berlioz 

cla-ta , = 0.3, kg ^ ’^ 2 ’' ^“8®® 

to 1.4- . Therefore , Sj^(W) is approximately Inyersely 
proportional to the particulate mass only when the latter is 
relatively ejuite large (i.e. , 102 '* ’'h ‘ ^ould he 

mentioned tlu-t the present method for deteimlnlng the 
dependoneo of on W is highly convenient , in that the 

more direct techniques for evaluating the data on the rate 
par® et era as a fuiiotlon of hold ^ of the particulate mass 

• j -hi tr iyT*fhflter exDerimental and computatio- 
would en-tail considerably greater e^jpeixm 

nol efforts - 


10.2 iiffGct of Ball Size 

+ •^/^n nf the nublished data of Kelsall 
Careful inspection of tne p 

. 4 .n ciweests that there exists a unique 

et al, I 94 I ’ suggests 




Fig-10.3 Cunnulativ <2 fraction retained . on the stated size as a function of t 
scaled grinding time for six different values of the particulate mass 
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relationship between the rate parameter Sj_(B) and the ball 
size B ill the following manner 


Si(B) 


Y f S -m . . -j 


( 10 . 11 ) 


where is the maximum value of the rate paraneter for 

the i— th size interval corresponding to the optimal ball size 

* ‘-,0 . 

, and u is a constant inteipreted as the smallest ball 
size ncceccar.y for any meaningCul grinding of the smallest 
particle sine of interest . Ihe function T is a unimodal 
frequency f laiction . In order to obtain a general relationship 
it is iioccsGary to express S and B as functions of seme 


characteristic size of the particles lying in a siae interval i. 
When sieve sizes differ by a constant factor 6 , we are free 


to ret>rorunit this characteristic size by the lower sieve size. 


upper sieve size , arithmetic or geometric mean of the lower 
and upper sieve sizes ♦ In each case the ratio of the 
characturir^tic sizes for adjacent size intervals is the same , 
Kie characteristic size chosen here is the size of the 
openings of the lower sieve , denoted by x^^ for the i-th 
size intei*val . 

c* 

Let the functions describing the variations of 
and B* witis the sieve size interval be fomally represented 


as 

B* = 


(1C. 12} 
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and 


(10.13) 

0. (I0.ia) ana (10.,,) anac,. (10.11) 


Sj^(2i) = Ug (U^ (x^)^ 




] (10,14) 


vvhion i 3 the desired correlation lor the rate paraneter 
ac a f'aictjo.n of ball and particle size 

Por dtrionotration data lor quarts grou.nd in a 
oontieuo.ie ball mill wore read off from ilg . 15 of leisaii 
ot .ii . I parameters were reported as a function 

of five kill uises for 6 p.rtiole else internals , correspond- 
ing to iower oibo lliuits of 833 , 590, 418 , 295 , 208 and I 04 
mi oiraij a-.u ciove nizc ratio of Y2 . The hold-up of the 
particMlvtc liawterial in the ball mill varied with ball size 
as 20(0, 1'.02, 1390, 1237 and 1216 gm. for 0 . 5 , 0 . 75 , 1.0, 

1,.-. < atifi. 1 , ..lO iiio/i balln , respectively . Assuming that the 
rate j.'arn.; ot i.,rc v;erc inversely proportional to the hold-up 
of t:’!c b ^.rtioul..t 0 mass , the values rq)orted were 


nomciil.jul by heis.'ill et al, for a standard hold-ip mass of 
1350 gjn, loi view of the conclusions arrived at in the 
previous ;,:cction , it is likely that the corrected values 
of the paraioters for 0.5 inch balls may be in appreciable 
error , if the relative magnitudes of the terns k^ and kgW in 
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Eq, (10.10) were of the satne order , However , for 0,75, 1.0, 
1.25 and 1,5 inch balls , the maximum deviation in the hold-up 
mass a.ooi.rl the standard hold-up of 1350 gn. was about 11 per- 
cent only , Therefore , the normalizing factor W/1350 
(ff = 1502, 1390, 1237, 1216 gn. ) applied by Kelsall et al, 
should rive a]iproxlmately the same results as those obtained 

ki + kg ¥ 

by Lijiny tinj factor s as per the correlation 

in Eq. (10,10)* . 

TliC plots of the data of Kelsall et al. in ELgure 10.4 
sho\; that Eqs. (10,12) and (10,13) can be written as 


•b. 

1 

* 


^'l ■*“ ^2 ^i 




K, X, 
3 1 


(10.15) 

(10.16) 


where K-, xin and K-. are constants which can be readily 
computed fiV-m; the curve . Eigure 10.5 shows the graphical 
forii of t.u. function ? in Eq. (10,11) . On inspection it 
turno cut that the curve can be fitted to the Gaucliy 

dintributioi. function | ^24 I Ps ’ 

K. 


V(“) 


( 10 . 17 ) 


+ (log P)' 


* Ir. oTtc of Uerlioz data the two oorreotlon procedures would 
reoul.t ih r>axh.um 5 poivoert aifferehoe In the values of 
tlia rate paraiieters over the range 1200 to 1500 ga.hold 



OPTIMAL BALL DiA , 3 , inch 
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Fig, 10. 4 Optimal bail diameter and the product of 

maximum value of rate parameter and optimal 
ball diameter as functions of particle size- 




'"o 1 i ”' 3 V" ' 5' 6 ■ 

ib-b°)/(b''-b°) ■ ' . , '■ 

Fig 10. 5 The ratio of rate parameter to its maximum 
value as function of variable p showing the 
self- f'rcscrving curve. ^ 
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where slid are parameters of the distribution and 




P - “i — 

- B° 


(10.18) 


Since , for 
and (1C. 17) 


-8 , S = S and p = 1 , hence , from Eqs. (10.11) 


K 5 - (10.19) 

Oonjvj.’ u I'.'y? '-'S yhown in figure 10.6, a plot of S /S as a 
function 01 (log p) is a straight line with slope equal to 
■|/K rind iitercept 1.0 . The final expression in Eq. (10.14) 

/ 4 


j]Q\V jej C 


Si (‘0 




( 10 . 20 ) 


[^1 + ^2 P4 ^ 

Figure 10.7 the family of curves has been back 
oaLcuj ot,(;d from Eq. ( 10 . 20 ) und compared with the experimental 
data , With; the exception of the data points pertaining to 
O.b iacl i-Jio , vvhich was to be expected for reasons mentioned 
above , tnc at:ruaaout in general is quite satisfactory. 


10. 3 Ciu'sc I ons 

7 /<,- iittve shown, at least for the grinding sys 


invuoti gated » 
paranet c.r!.’.! for 


that to a first approximation the rate 
different size intervals are influenced by the 



I', i . -,)(/ r. « rnn..ton 


O H J3 


Inverse of the ratio of rate parameter to its. 
maximum value as a function of (log p)^ for 
estimation of constant K4 in the cauchy distri 





Ui. / 1 I LJ 

0.25 0.5 0.75 1.0 1.25 1.5 

BALL DIAMETER, inch 


Fig -10 .7 Comparison of experimental and predicted 
values of rate parameters for different 
ball sizes- 
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particulate niass to the aane extent, in the manner of 

Bl. (10.1) . The parameter 3^(1?) Is related to the partleulate 

mass W ae Ih Bq. (10. 10) , which is found to he valid over 
both the unsaturated and saturated regimes , eioept when the 
mill is hi silly starved . The inverse proportionality between 


the r..te p..raaeter and particulate mass , as suggested by 
earlier au^hora , is valid only in the limit for hi^y 
saturated ;?:illj? , D«ie to lack of ^propriate data , it has 
not eec.. pciicib: c to determine if the function H(W) depends 
on tile J,J size or not , although the latter was Implicitly 


ass'-fcicd f. r tiie .rinalysis of the data of Kelsd.1 et al. 


T!u. sdLf preserving curve in Mgure 10.5 , where 
s/s* is I'lottiiri as m f<uiction of (B-BP)/'B*-JB®),..‘2;3implifies 
conoid 3\d)iy oin: sonroi: for a reasonably accurate correlation 
for the .•-•ri'ect of the ball size . Ihe choice of the Oaiichy 
distrih!i'M.L. i. if! oor.vciiJ.ent for gr^hical computation of its 
par^r.cucr , but it is by no means unique « It is possible 
that othiir tii retributions can be fitted to the seif-p reserving 
curve • 


r'inally, it is recognised that the nature of the 
fiuiction.'e H(iy) , , fJg and V are specific to the grinding 

systai.-: invvstitjritcd and there is no 'a priori’ reason. why 
these shotaa rtwain invariant for other systans' .. .Wever ., 
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WG ivive ud details for checking the assumptions made 

in Bqs. (10,1) , (10,2) and (10.11) , and given simple 
gr* ohical -rccoduics for determining the enpirical correlations 
for tlic cif^otG of p 'I'ticul ate mass and ball size on the rate 
par'miv I", rr-. fro::; suitaulo experimental data . 



CHAPTER - 11 


SClvIMARI AKD CONGLUSIONS 

A new solution has teen obtained to the fundamental 
integro-diff erential equation of batch grinding for selection 
and breakage functions of the form given in Bqs. (4.1) and 
(4.2) , rei^ectively. Ihe solution given in Eq. (4.13) is the 
most general closed form analytical solution available at 
present. It describes the complete trajectory of particle 
spectra in time, ^art from a small initial period. Unlike 
some commonly used s el f-p reserving size distribution equations 
such as Gate-Gaudin-Schumann distribution and Rosin-Rammler- 
Bennett distribution \ 13 \ , Eq. (4.15) has anon-self-preserving 
structure • 'JOhis point has been illustrated by conputer 
simulation experiments for two different forms of the basic 
function g(x) in Bqs. (4.1) and (4.2) . 

2, Using a formal mathematical discretization procedure, 

exact expressions have been obtained for the parameters of the 
two size-discretized models in Eqs. (2.86) and (5.10) in terms 
of the three basic functions S(x), B(v,x) and M(x,t) of the 
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size continuous, linear and time - invariant mathematical 
models Equations (5.7) , (5.8) , (5,15) , (5,15) and (5.1?) 
clearly show that all the size discretized parameters are 
uniquely interrelated through these functions and therefore, 
cannot be assigned any arbitrary values. Moreover, because 
of the dependence of these parameters on the mass density 
function M(x,t) , those parameters are inplicit functions of 
grinding time and hence are non-unique, 

The product term S(v) B(v,x) occurs r^eatedly in 
the analysis of discrete size models. It was felt that it 
deserves a na;ae of its own. In this work it has been called 

as 'Selekage Function* and denoted by L(v,x) . Its important 
properties have also been delineated, 

5. It has been siiown that by specializing the 

functional forms forthe functions M(x,t) and L(v,x), some 
simplified size discretized models can be obtained. For 
sufficiently narrow sieve size intervals, it is a fair 
^proximation to assume that the particle size distribution is 
uniformly distributed in a piece-wise manner, over the 
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individual discrete size interval size r^g%s. Ihis gives us 
the Bass niodel, \'/itn time invariant discrete size parameters. 
With examples it has been shown that even when the sieve 
size ratio 6 is as small as '^Y'^ , the inner breakage parameter 
is not iiegLigLbl^y smctllj as has been assumed by some oi the 
investigators It the selekage function does not exhibit any 
dependence on the parent particle size v i.e. I(v,x) = L2(x), 
the RSi' model io obtaiiiod. This type of selekage function 
generates ROvSiii-Ban.ilcr type or Rosin-B.ammler distributions. 
The rate and brealc; ge parameters are time - invariant while 
the selection and distribution parameters are time dependent, 
For illustration, coir.putod results on time variation of 
b- .(t) for Lo(x) “ Ax have been graphically displayed for 

X X 

some ooluotc.d dis.oroto size intervals. The results show that 
at Icjint for uoor.v; si. so int erv -al s, the approximation of time 
invnri;:ait iuijc.r hruakage jiurameter is not v.alid. When the 
selekage tbsictihj is seperable i.-o, , l(Vi,x} = I>[(v) L2(x) , 
the SS? model is b t ;;..i n cd . In this model B . is time- 

Xy J 

invariant aiid rest of the parameters are time dependent. Time 

variation of b. .(t) has been coiiputed using 

1 1, r ct-p p 

Kapur's similarity solution for I(v,x) = A v x .It has 
beem concluded tuat the extent of time variation is 
determined by size interval width as well as the dependence 
'^f the selekage function on the parent particle size v . If 
the selekage function is weakly dependent on v , wider size 
intervals can bo eniployed in the analysis, with time invarian 
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parameters as a reasonable approximation. It has been further 
shown that when the selekage function I(r,x.) is a monotonically 
increasing function in size v , the disappearance kinetics plot 
for the discrete size interval i tends to concave upwards, and 
when it is a monotonically decreasing function in v , it tends 
to concave downwards v/ith increasing grinding time 

4. dt has bocij shown that, in general, it is not possible 

to realize aii order kinetics concept for a discrete size 
interval. Only in case of J3ass and RSI’ models? the dis^pearance 
kinetics of a discrete size interval is first-order, in the 
conventional sense. The phenomena of zero-order production 
of fines and Rosin-Rammler type grinding kinetics, are 
coup at i able witii only RSP model. The Alyavdin distribution 
need not nccoss<.irily arise from any inherent time-d^endence 
of the selection ftjncticn in Bq. (2.65) as assumed earlier. 

It has been shown , using synthetic size distribution data, 
that the Alyavdin equation could be a reasonably close 
spproximatc solution to theintegro-diff erential equation of 
batch grinding for a class of tima- in variant selekage 
functions . The brook age parameters difference- 

similar (norraaiisiiblo) in the Bass model if and only if S(v) 
is proportional to v^ , B(v,x) can be represented as B(x/v| 
and the sieve size ratio 6 is constant. The breakage 
parameters fjre difference - similar in the RSB and SSP models* 
also^ provided L(v,x) is proportional to . The rate parameter 
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is proportional to x“ in Bass and Ssj models only, pro 
S{t) Is proportionsl to v“ aai B(v,x) Is of the form B(x/ 

The phenomena of difformoe - siiailar breakage parametoR 
a power law for size variation of the rate paraaeters are 
always coupled. Tliosc results have clearly sstablished tha 
a giver) discrete size model possesses concurrent properties 
none of wirdch inunt be violated, otherwise serious inner 
contradictierji'. will ocvclop. 

5. A novol approacli has been presented for identificat 

of a linciox <md time im-nriojirt grinding systan in which size 
reducti .ui occur.s by diippiiig, impact fracture as well as 
abraci ur, A :rjill r tiictici) 0(v,x.<vt) has been defined which 
represents t!;u uizo distribution of the progeny particles 
when a unit impuicc of size v is ground for a fixed time 
interval t , Starting-; fri.'m a generalized polynomial in x 
and V , a suitable uyq'jreosion has been obtained for the mill 
function whieii is e .i.rsiatent with the pljysical constraints 
on the mill fiinctivoi, The constants in this expression have 
been estimated uojii/*- tiu, standard least squares method of 
fitting expurimentri! dr.tca. Two sets of batch ball mill 
grinding data f jr d, IcKutc have been employed for illustration 
of the ideaitif ication scheme. In both the cases, the fit 
to the c*xperiiiient.:d data has been found to be very good. The 
role of the vcucioua grinding mechanisms has been discussed 
in light of the nature of the estimated mill functions, and 
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been discussed . 


An eXcicfe olgoi’i'thiii h.as been developed for back 
calculation of tlio rate and breakage parameters which makes 
^ 3 e of Reid’s s ;lution to the size discretized batch grinding 
equation a ’0-1' vrriehle suppression technique . The 

grinding data required is single size feeds ground for a 
fixed tine inti.U'val t , which can also be generated with the 
help of ostii.iated i.iill function . Using this method it has 
been possible to cn:puto the time variation of the rate and 
breakage paraxxterc due tu their dependence on the size 
distribution ;f th.; particles within the individual discrete 
sise intervals. VVith a view to providing a standard reference 
for c mp-irisan jf dit'fcrciit grinding systons, anew set of 
par.ar..utcrn - 'reforonce parameters' have been defined. A 
coi%>arison f the caputed reference parameters has been made 
with the r;ppr, xiuiatc c-stirnates of rate and breakage parameters 
repozted in the literature fwi' the same grinding systems. The 
=h.,'.vn i.i prinoiplo the mill funotion ^proaoh 

and the pr,;p.iK.;d taaitific-ition scheme Is quite sound. 


7. Systomr.-io rruphical pipcedures for deteimlnaticn of the 

• ■ t .'....I ‘ in*!'’ for the variations in time independent 

empirical no rruiauitwii. xor tue 

. . / • -•odel) with particulate mass and the size 

rata parx.atarn ( ouci; .v 

, i.;.. i...,’; ■ f* have been developed and successfully 

of tlio j:rin4iriy ;:u.nlia satib nuvo u 

demonr.triitf?!! fur tuo acfcuiO. grinding oyst 
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APPENDIX I 


A Relationship Between Characteristic Particle Size and the 
Pirst Moment of Particle Size Distribution | 

K^ur j 14 I has shown that under the restrictions given 
in Eqs. (2.56) and (2.57) a similarity sc-lution to the 
integro-diff erential equation of batch grinding is 

= 'hp) z(|) (1.1) 

where 

^ ( 1 . 2 ) 

The k-th moment of the particle size distribution at time t 
is given by 

= / M(x,t) dx (1.5) 

Thus, combining Eqs. (1.1) , ( 1 , 2 ) and (1.3) gives 


11 

J [ / 

( 1 . 4 ) 

or 

, k 


It 


(1.5) 

where ^ is 

a dimcsnsionless constant . Prom Eqs. 

(2.4a) and 

(1.5) it follows that 




(1.6) 


Hence, it is proved that the characteristic particle size x is 
proportional t.i the first raomc-nt of particle size distribution. 
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(11*3) 


In principle, this solution sliould be equivalent to the 
one obtained directly from Eq. (2.59) when a is replaced by ad 
and ^ is replaced by pd. Substituting for in terms of 

j) 

from Eq. (11,4) > it can be cc^nfirmed thfit the solution 
in Eq. (11,5) really reduces to the self-similar form , and 
is identical to the one obtained directly from Eqs. ( 2 , 59)5 
(2, 60) and (2. 61 ) . 
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i^PPEUDII III 

F1...W CSi.-ort £■ r Bo.ok Calculation of the Parameters in the 
Sise jOLscrotined 
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